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' , Abstract: We present brief synopses of supersymmetric models where either the neu- 

QJ I traUno composition or its mass is adjusted so that thermal relic neutralinos from the Big 

' Bang saturate the measured abundance of cold dark matter in the universe. We first review 

minimal supergravity (mSUGRA), and then examine its various one-parameter extensions 
^ ' where we relax the assumed universality of the soft supersymmetry breaking parameters. 

. Our goal is to correlate relic-density-allowed parameter choices with expected phenomena 

^T) . in direct, indirect and collider dark matter search experiments. For every non-universal 

I model, we first provide plots to facilitate the selection of "dark-matter allowed" parameter 

I space points, and then present salient features of each model with respect to searches at 

QQ ' Tevatron, LHC and ILC and also direct and indirect dark matter searches. We present 

■ benchmark scenarios that allow one to compare and contrast the non-universal models with 

one another and with the paradigm mSUGRA framework. We show that many implications 
^ . about sparticle properties and collider signals drawn from the analysis of the relic density 

^ . constraint within mSUGRA do not carry over to simple one-parameter extensions of the 

mSUGRA framework. We find that in many relic-density-consistent models, there is one 
(or more) detectable edge in the invariant mass distribution of same-flavour, opposite sign 
dileptons in SUSY cascade decay events at the LHC. Finally, we scan the parameter space 
of these various models, requiring consistency with the LEP2 constraint on the chargino 
mass, and with the observed relic density, and examine prospects for direct and indirect 
dark matter detection. We find that in a large number of cases the mechanism that causes 
the early universe neutralino annihilation rate to be large (so as to produce the measured 
relic density) also enhances the direct detection rate, and often also the rates for indirect 
detection of neutralino dark matter. 



Keywords: pupersymmetry phenomenology, Neutralino dark matter, Collider and dark 



matter signals 



Contents 



|l]. Introduction and framework 



I 



[2]. Brief synopses of SUSY models with neutralino dark matter 



i 
i 



2.2 



The mSUGRA model 

Models with scalar mass non-universality 

2.2.1 Generational non-universality: normal scalar mass hierarchy 

2.2.2 Non- universal Higgs mass: one extra parameter case 
Non- universal Higgs mass: two extra parameters case 



2.2.3 



2.3 



Models with non-universal gaugino masses 
2.3.1 Mixed wino dark matter 



2.3.2 Bino-wino co- annihilation (BWCA) 

2.3.3 Low IM3I dark matter: compressed SUSY 



2.3.4 High IM2I dark matter: left-right split SUSY 



Illustrative benchmark cases 

General characteristics of relic-density-consistent models 

!l| Implications for collider searches 
2| Implications for (g — 2)^ and BF{b — > 57) 

Implications for direct detection of dark matter 

Implications for indirect detection of dark matter: neutrino telescopes 
Implications for indirect detection of dark matter from halo annihilations 








11 

14 
16 
18 
19 

m 



27 
30 

m 



m 

40 
42 



Isl Summary and concluding remarks 



46 



1. Introduction and framework 

An abundance of evidence arising from a variety of cosmological measurements shows that 
most of the matter in the Universe is not baryonic, but rather composed of massive neutral 
stable (or at least extremely long-lived), weakly (or super- weakly) interacting particles. 
Since none of the particles of the Standard Model (SM) have these properties, the existence 
of this so-called dark matter (DM) in the universe provides unequivocal evidence for physics 
beyond the SM. 

Cosmological measurements severely constrain the abundance of DM: combining the re- 
sults from the WMAP Collaboration with those from the Sloan Digital Sky Survey gives |T[ 

J^DM/i' = 0.11110;°^^ (2^7), (1.1) 

where ri = p/pc with pc the closure density of the Universe, and h is the scaled Hubble 
parameter, h = 0.73 it 0.04. While the mass density of DM is rather precisely known, 
the identity of the DM particle remains a mystery. One class of candidates - thermally 
produced weakly interacting massive particles (WIMPS) - are especially appealing in that 
they naturally occur in a variety of well-motivated models, and further, because their 
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relic density today (which can be rehably computed) is found to automatically have about 
the observed magnitude, provided the WIMP mass is of order the weak scale: mwiMP ~ 
100 GeV. Of course, we should always bear in mind that, like visible matter, DM may 
consist of several components, so that, in standard Big Bang cosmology, (1.1) really implies 
an upper bound on the density of any single component. 

Softly broken supersymmetry (SUSY), with a SUSY breaking scale below 1-2 TeV, is 
highly motivated for a variety of theoretical as well as experimental reasons |Q, Q. SUSY 
models with conserved i?-parity include a stable, massive weakly interacting particle - 
the lightest neutralino Z\ in many models - which is perhaps the prototypical thermal 
WIMP Q. In any supersymmetric model with a stable neutralino, the neutralino relic 
abundance can be reliably calculated as a function of model parameter space [^. The 
result depends inversely on the thermally averaged neutralino-neutralino annihilation and 
co-annihilation cross sections, integrated over time from freeze-out to the present day. Once 
the parameters of the model under study are known to match the measured relic abundance 
( |1.1| ) 1^], then these select parameter space regions can be checked for phenomenological 
constraints from low energy measurements and from non-observation of new physics signals 
in the LEP and Fermilab Tevatron data. Implications for the on-going Tevatron run as well 
as for experiments soon-to-begin at the CERN LHC, and possibly at a TeV linear electron- 
positron collider in the future can be examined. Likewise, predictions can be made for 
rates of direct detection of relic neutralinos via scattering on nuclear targets, or rates for 
indirect neutralino detection, either via signals from neutralino annihilation in the solar 
core, or via galactic halo annihilations which can give rise to gamma ray or anti-matter (p, 
or Z)) signals. 

Most analyses of neutralino dark matter have been carried out in the context of the 
minimal supergravity model - mSUGRA Q where SUSY breaking, which occurs in a 
hidden sector, is communicated to the observable sector via gravitational interactions. 
The universality of soft SUSY breaking (SSB) parameters, renormalized at a scale Q ~ 
Mq,\}'y — Mp is the hallmark of this framework. Specifically, one assumes that the mediation 
mechanism induces a common mass parameter mg for all MSSM scalars, a common gaugino 
mass mi/2 foi' gauginos, a common trilinear SSB parameter together with a bilinear 
Higgs scalar mass 6, in the effective MSSM Lagrangian, with parameters renormalized 
at Q = Mgut. It is also assumed that the dimensionful SSB parameters all have the 
magnitude of the weak scale. The large top quark Yukawa coupling drives the celebrated 
radiative electroweak symmetry breaking (REWSB) mechanism, and automatically leads 
to the S'[/(3)c X C/(1)em symmetric vacuum over a significant portion (but not all) of the 
model parameter space. The GUT scale SSB parameter h can be traded for tan/3, the 
ratio of Higgs field vevs, while the magnitude (but not the sign) of the superpotential mass 
parameter [i is fixed by the observed value of Mz- The mSUGRA model is thus completely 
specified by the parameter set 

mo, mi/2, ^0, tan/3, sign{fi), (1.2) 

along with the value of the top quark mass nit- Except where explicitly mentioned, we 
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fix mt = 171.4 GeV, in accord with recent top mass measurements at the Fermilab Teva- 
tron |]. 

Unless sparticles are very light (~ 100 GeV) the generic value of the relic density in 
mSUGRA (as well as in many other SUSY models) tends to be well in excess of the observed 
CDM relic density ( |1.1D . As a result, only special regions of the mSUGRA parameter space 
where the annihilation rate for neutralinos is enhanced are compatible with the measured 
value of the relic density.^ In early work on the mSUGRA model, the low mo, low mi/2 
region (so-called "bulk region"), where sparticles are indeed very light, was favored Q in 
that neutralino annihilation into leptons via light t— channel slepton exchange occurred at 
large rates, leading to relic densities f^^i^^ ~ 0.3—1. The rather lower measured abundance 
in ( |1.1| ), however, favors even lower values of mo and mi/2) resulting in considerable tension 
with the negative search results from LEP2 for chargino and slepton pair production. 
Within the mSUGRA framework, the remaining relic-density-allowed regions consist of: 

• The stau-co-annihilation region at low mo and low-to-moderate values of m^/2 where 
uif-^ ~ m^^ , so that neutralinos can co-annihilate with staus ||l^ in the early universe. 



The hyperbolic branch/focus point (HB/FP) region |11] at very large mo values where 
becomes small so that Zi becomes a mixed bino-higgsino state. In this case, ZiZi 
annihilation to WW, ZZ and Zh via the Zi higgsino component is enhanced in the 
early universe. 

The Higgs funnel region at large tan/3 ~ 50 [12|, where 2m^^ ~ niA, so that neu- 
tralinos can annihilate at an enhanced rate through the (wide at large tan /5) A (or 
H) resonance. An /i-resonance annihilation strip can also occur at low m^/2) where 
2m^^ ~ m/i Q . 

The top squark co-annihilation region at large negative ^o values where m^^ ~ m^^ 



so that Zi can co-annihilate with ti particles |14]. 



The regions of mSUGRA parameter space leading to a neutralino relic density in agree- 
ment with ( |1.1| ) are all near the edges of theoretically (or in the case of h resonance annihila- 
tion, experimentally) allowed parameter space, which can lead one to question whether the 
mSUGRA model might be disfavored by the measured neutralino relic abundance. In this 
vein, many authors have examined SUGRA-type models but with non-universal soft term 
boundary conditions at Q ~ Mgut- R is appropriate to note here that unfettered non- 
universality of soft terms generically leads to the occurrence of flavor changing processes 



at levels far beyond experimental limits [15|. With this in mind, we work in a simplified 



parameter space wherein there exists degeneracy or near degeneracy of first and second 
generation scalar masses, leading to a suppression of FCNC processes via the super-GIM 
mechanism. At the same time, in order to maintain the obvious success of gauge coupling 
unification, we must assume that the correct effective theory between the weak and GUT 



^In our analysis, we are assuming that thermally produced neutralinos in the standard Big Bang cos- 
mology make up the DM. While it is possible to get around these assumptions, we feel that an examination 
of the conceptually simplest scenario that does not invoke additional hypotheses warrants special attention. 
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scales is the MSSM, or the MSSM augmented by gauge singlets, or the MSSM together 
with additional matter in complete multiplets of SU{5). 

With these considerations in mind, we will assume that: 

1. In the interests of minimality, while maintaining the successful predictions of gauge 
coupling unification, that the MSSM is the correct effective field theory between 
Mweak and Mqut- 

2. The REWSB mechanism leads to an SU{3)c x C/(1)em symmetric ground state; i.e. 
electric charge and color gauge symmetries are not spontaneously broken. 

3. CP violating phases in the SSB parameters are suppressed so that supersymmetric 



contributions to C-P violating processes are sufficiently small [16| 



4. There is a near-degeneracy of SSB of the first and second generation sfermions so that 
SUSY contributions to flavor-violating processes is automatically suppressed. We 
do allow some non-degeneracy between third generation scalars and first or second 
generation scalars. 

5. i?-parity is conserved so that the lightest supersymmetric particle (LSP) is stable. 

6. The gravitino - which in models with gravity-mediated SUSY breaking naturally has 
a mass of order M^^ak ~ is not the LSP, which we take instead to be the lightest 
neutralino. For a discussion of the possibility that a gravitino LSP is the DM, see 
Ref. 11]. 

In the spirit of our earlier discussion we need to relax the theoretically least well- 
motivated universality assumption that underlies the mSUGRA framework in a controlled 
manner (to avoid large flavor-violating couplings) and explore non-minimal SUGRA models 
with an expanded parameter space. We could, for instance, consider a non-minimal model 
where we independently vary. 

Ml, M2, M3 (gaugino masses), (1.3) 
mo(l,2) (common first/second generation SSB matter scalar masses), (1-4) 

mo(3) (common third generation SSB matter scalar masses), (1-5) 

n^'n^, "^//^ (non — universal SSB Higgs mass parameter), (1-6) 

At, Af,, A-r (non — universal third gen. ^ terms), (1-7) 

tan/3, (1.8) 

sign{fj,). (1.9) 

A different but equally reasonable option may be to require common masses for matter 
scalars with the same gauge quantum numbers but allow intra-generation splittings. In this 
case, we would have common mass parameters rrig, mfj, mj^, m\ and m?^ at Q = Mgut- 
In the extreme case, the matter scalar masses can be further broken down into specific soft 
term masses toq., 'rrifj,, "rnfj,, '^iji ^1;^ where i = 1 — 3 for each generation. These options 
have been explored elsewhere [18, IS], and will not be discussed further in this paper. 
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One way to proceed is to perform scans over the much larger non-minimal SUGRA 
parameter space and search for solutions which satisfy dark matter (and also other) con- 
straints. While this approach has the virtue of being unbiased in the scanning, it is practi- 
cally difficult to implement. Moreover, when a large number of free parameters are varied 
simultaneously, it is frequently difficult to draw insights into the associated dark matter and 
collider phenomenology that follow. Instead, many groups [20, 21, 22, 23, 24, 25, 26, 27| 



have examined the impact of relaxing the underlying universality of the mSUGRA model 
in a controlled way, by allowing non-universal parameters only in one sector of parameter 
space at a time. For instance, we may consider the mSUGRA parameter space augmented 
to accommodate non- universal gaugino mass parameters [^] as in ( |1.3D , or instead ex- 
tended to allow Higgs boson SSB mass parameters to be different from matter scalar mass 
parameters as in ( |1.6| ), but not both. Other directions in the parameter space of the non- 
minimal SUGRA models can be similarly explored. These extensions generally require 
augmenting the mSUGRA space by just one (sometimes, two) additional parameter that 
is adjusted to yield agreement with the observed DM relic density. The phenomenological 
implications of the extended model as a function of the remaining mSUGRA parameters 
can be readily examined, and directly compared with the paradigm mSUGRA framework. 
This approach has led to new insights and to exciting new possibilities for collider and 
dark matter phenomenology that can be expected in models with non-universal soft SUSY 
breaking terms. 

Examination of these simple one-parameter extensions of mSUGRA leads to another 
important pay-off. Since, as discussed above, analyses of the relic density constraint in 
mSUGRA force parameters to be in the bulk region, the stop or stau co-annihilation 
region, the Higgs funnel region or the HB/FP region of parameter space, many groups 
have inferred that at least one of the following must hold: 

1. Sfermions have masses ~ 100 GeV (bulk region), and so must be accessible at the 
LHC. 

2. There is at least one charged sparticle close in mass to the LSP, so that this should be 
accessible at the LHC, unless the LSP is so heavy that the hard- won gauge hierarchy 
is again destabilized (co-annihilation). 

3. The additional Higgs scalars of the MSSM are relatively light with rriA ~ 2m^^ 
so that these can be searched for at the LHC, which requires large values of tan (3 
where sparticles preferentially decay to third generation quarks and leptons (Higgs 
resonance region). 

4. The lightest neutralino has a significant higgsino component, which is possible only 
if mo is so large that squarks and sleptons are essentially inaccessible at the LHC 
(HB/FP region). 

It is imperative, of course, to check the robustness of these "predictions" to minor variations 
of the assumptions underlying mSUGRA before drawing broad conclusions about what 
the relic density determination implies for experiments at the LHC, as well as for direct 
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and indirect detection of DM. Our study of the various extensions of mSUGRA naturally 
permits this. 

In this paper, we have two broad goals. The first, in Sec. ^, is to present an overview 
of a number of different models, wherein by tuning one additional parameter beyond those 
of the mSUGRA model we can match the predicted neutralino relic abundance with 
Models wherein the composition of the neutralino is adjusted to obtain the measured relic 
abundance are referred to as "well-tempered neutralino" models (WTN) [^. We also 
examine several models wherein neutralino or other sparticle masses are adjusted to obtain 
the correct relic abundance of dark matter. 



In each case, we present i.) motivation, ii.) the parameter space, and selected param- 
eter values that allow the reader to generate spectra and collider events for the particular 
model. We also comment on the salient features of Hi.) collider and iv.) dark matter 
search phenomenology associated with each particular model. For our analysis, we adopt 
the SUSY spectrum generator Isasugra, a part of the event generator ISAJET 7.76 pO|. 



For any given parameter set satisfying the DM relic density constraint (1.1), the sparticle 
mass spectrum and associated neutralino relic density and direct and indirect detection 
rates may be calculated, and associated collider events may be generated for the Tevatron, 
LHC or ILC colliders. To facilitate comparison with the paradigm mSUGRA case, we first 
present an updated overview of allowed regions within mSUGRA. In Sec. ^, we present 
some benchmark cases where the spectra and some results from these various models are 
explicitly compared with the corresponding situation in the mSUGRA case. 

Our second goal, presented in Sec. is to extract several general results from scans 
over the models examined in Sec. to gain an idea of some of the features relevant to 
collider and dark matter searches that might be shared by many of these models. For 
instance, it has already been pointed out that the subset of these models which resolve the 
dark matter relic density problem via mixed gaugino/higgsino dark matter (i.e. models 
with a WTN) collectively have neutralino-nucleon direct detection scattering rates that 
asymptote around ~ 10~^ pb|^l], putting them within reach of direct dark matter search 
experiments currently being mounted, such as SuperGDMS, LUX, Xenon-100, WARP and 
mini-CLEAN. Also, models with non-universality where the composition is tempered to 
yield the observed relic density, or where agreement with ( |1.1D is obtained via bino-wino 
co-annihilation, tend to have a neutralino mass gap — m^^ smaller than Mz, so 
that three body decay modes dominate the Z2 branching fraction. Unless the leptonic 



branching fraction of the neutralino happens to be strongly suppressed |32|, this then 
yields an observable mass edge in the dilepton mass spectrum at m{i^i^) = m^^ — tti^^, 
which serves as a good starting point for sparticle mass reconstruction in gluino and squark 
cascade decay events at the CERN LHC [33|. 

We conclude in Sec. I with a summary of our results together with some general 
comments. 
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2. Brief synopses of SUSY models with neutralino dark matter 



2.1 The mSUGRA model 

We begin by presenting updated results on the allowed parameter space of the minimal 
supergravity model. The mSUGRA model is completely specified by the parameter set 



O). To calculate the sparticle mass spectrum, we use ISAJET 7.76 |30|. The relic 



density is evaluated via the IsaReD program [34|, which is part of the IsaTools package. 



IsaReD evaluates all 2 — > 2 tree-level neutralino annihilation and co-annihilation processes 
and implements relativistic thermal averaging in the relic density calculation. 

For our first results, we show in Fig. || the mo vs. mi/2 plane for parameters = 
0, tan/3 = 10 and > with a) nit = 170 GeV, b) nit = 171.4 GeV and c) nit = 
175 GeV. The red-shaded regions on the left are excluded because fi becomes the LSP, 
while the red-shaded regions on the lower right are excluded due to a failure to meet the 
EWSB minimization conditions. The blue-shaded region is theoretically allowed, but is 
experimentally excluded by LEP2 searches for chargino pair production where we require 



> 103.5 GeV [Q. The negative results of Higgs boson searches at LEP2 |36] require 
that the SM Higgs boson is heavier than 114.1 GeV. This limit can be translated to a 
lower limit on the MSSM Higgs boson mass. While h ~ -ffsM if innA is large, in general the 
bound on rrih depends on MSSM parameters, including CP violating phases that we have 
ignored in our analysis. The evaluation of m/j is also uncertain to about ~ 3 GeV due to 



missing two-loop corrections |37]. For these reasons, we do not include any bound on ruh in 
the LEP2-excluded blue region, but only show the boundary of the region m/j < 110 GeV 
by the magenta contour (lower-left) in the figure. The green regions have a neutralino 
relic density in accord with (^]^): 0.094 < ^^_^h'^ < 0.129. In the yellow regions, however, 
^Zi^^ < 0.094, so that an additional component of dark matter particles is necessary to 
saturate the observed DM relic density. The remaining unshaded regions all have h'^ > 
0.129, i.e., they give rise to too much dark matter: thus, these are excluded in standard 
Big Bang cosmology. We also show contours of gluino and first generation squark masses; 
these contours hardly change under variation of Aq, tan/3 and sign{jjL), except at the level 
of one loop corrections and L'-term contributions to their masses. 

The hard-to-see green/yellow region adjacent to the fi-Zi region shows up as a very 
narrow sliver where the fi — Zi mass gap is small enough so that stau co-annihilation 
occurs at a large rate. This region in fact appears jagged only due to the resolution of 
our parameter space scans. One can also see the HB/FP region - where |//| becomes 
comparable to the SU{2) and U{1) gaugino masses and the Zi becomes mixed higgsino 
dark matter - adjacent to the EWSB forbidden region as the wider green/yellow shaded 
region at large mo, starting at mx/2 ~ 300 GeV, which corresponds to the turn-on point 
for ZiZi W^W~; for lower mi/2 values, this annihilation channel is closed, and the 
neutralino annihilation rate (via Z* exchange) generally becomes too small to bring the 
relic density into accord with (|l.l| ). 

While these three frames for the different m^ values are qualitatively similar, the main 
effect of mt variation shows up in the location of the EWSB excluded region, and hence 
the location of the adjacent HB/FP region: on the low side of the allowed nit range, the 
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mSUGRA : tanp-10, \=0, \l>f) 



a.) ra = 1 70 GeV b) m = L 7 L .4- Ge V 




c) m=n5 GeV 




I S 3 -I S 

m„ (TsV) 



Figure 1: A plot of the mo vs. mi/2 plane in mSUGRA for Aq = 0, tan/3 = 10 with n > and 
a) mt — 170 GeV, b) rrit = 171.4 GeV and c) = 175 GeV. The red-shaded regions are excluded 
because electroweak symmetry is not correctly broken, or because the LSP is charged. Blue regions 
are excluded by direct SUSY searches at LEP2. Yellow and green shaded regions are WMAP- 
allowed, while white regions are excluded owing to ^^^^^ ^ 0.129. Also shown are gluino and first 
generation squark mass contours, as well as a magenta contour below which nih < 110 GeV. 



HB/FP region moves to mo values as low as 1.5 TeV, while at the high end of this range, 



the HB/FP region only starts when mo ^ 3 TeV [38 



In Fig. 1^, we show the mo vs. mi/2 plane for Aq = 0, fi > 0, nit = 171.4 GeV, and for 
six different values of tan /?. We see that for tan /? = 10 only the stau co-annihilation and 
HB/FP regions are DM-allowed. As tan/3 is increased, more and more parameter space 
comes into accord with ( |l.lD . Already for tan/? = 45, a small DM-allowed region appears 
at low mo and low mi^2- The reason is that as tan/3 grows, the b and r Yukawa couplings 
become large, causing niA to drop. Thus, ZiZi — > A* ^ bb, t^t^ becomes more and more 
important, even if one is not "right on the A resonance" [39|. The low mo, mi/2 allowed 
region grows even more at tan /3 = 50. At tan f3 = 52, the A- funnel annihilation region 
has just come into view on the left edge of parameter space. By tan /3 = 54, the A-funnel 
is extremely broad due to the large A width: F^ ~ 10 GeV (50 GeV) for low (high) mi/2- 
We notice at tan P = 55, a small red-shaded wedge invades the plot at low mo and mi/2- In 
this region, the value of m| < 0, signaling collapse of EWSB. For somewhat higher tan /? 
values, the entire parameter space collapses. 

In Fig. I, we show the same mo vs. m^/2 planes as in Fig. ^ for multiple tan (3 values, 
but this time for ^ < 0. This sign of fi is disfavored by the Muon g — 2 Collaboration 
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t, J.5 
1 



mSUGRA : Aq = 0, n > 0, = 171.4 GeV 

a)tanp^lO b) tai]p=30 




c) tanp=i5 



e) tajip^52 





d) tanp=50 




f) tan|3=55 




o<nti-< 0.094 
t 0.094 <m"<0.iZ9 
— m,= LLOGeV 



• Theoretically Excluded 

• LEP2 Excluded : ni~^ < i03.5 GeV 



Figure 2: A plot of the mo vs. mi/2 plane in mSUGRA for Aq = and various values of tan /3, with 
/i > and rrit = 171.4 GeV. The red-shaded regions are excluded because electroweak symmetry 
is not correctly broken, or because the LSP is charged. Blue regions are excluded by direct SUSY 
searches at LEP2. Yellow and green shaded regions are WMAP-allowed, while white regions are 
excluded owing toQ^ h'^ > 0.129. Below the magenta contour in each frame, m/, < 110 GeV. 



measurements of anomalous magnetic moment {g — 2)^ at low mo and low At 
high mo and mi/2 values, sparticle contributions to the muon QED vertex decouple, and 
the deviation from SM predictions is tiny for either sign of fi. We also see that the ^-funnel 
arises in the mo vs. rnxji plane at somewhat lower tan/3 values. The ^-funnel is actually 
narrower than in the /x > case, in part because the A width is narrower. We also see 
a bulge of incorrect EWSB beginning already at tan /? = 45, and growing so as to engulf 
nearly all parameter space by tan/3 = 55. 

We see from Fig. ^ that the Higgs funnel moves to larger values of mg as we increase 
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mSUGRA : Aq = 0, n < 0, = 171.4 GeV 

a)tanp^lO b) tEmp=30 




c) tanp=i5 





e) tanp^52 



f) tan|3=35 




0<nti-< 0.094 
t 0.094 <m"<0.iZ9 
— m,= LLOGeV 



• Theoretically Excluded 

• LEP2 Excluded : ni~^ < i03.5 GeV 



Figure 3: A plot of the mg vs. mi/2 plane in mSUGRA for Aq = and various values of tan/3, 
with /i < and rrit = 171.4 GeV. The red-shaded regions are excluded by lack of correct EWSB 
or by the presence of a charged LSP. Blue regions are excluded by direct SUSY searches at LEP2. 
Yellow and green shaded regions are WMAP-allowed, while white regions are excluded owing to 
^zj^^ > 0.129. Below the magenta contour in frames a)-c?) rrih < 110 GeV , while rrih > 110 GeV 
all over the LEP2 allowed region in frames e) and f) . 



tan p. To understand this, we take a point {mo,mi/2) in the funnel where ~ ^"^Zi ' ^^'^ 
examine what would happen if we increase tan /3 keeping the other parameters fixed. We 
first remark that because Zi is essentially a bino, ~ Mi remains essentially unaltered. 
The behaviour of ~ m'jj^ — rrij^^ is governed by how the evolution of the Higgs scalar 
SSB parameters is altered by the increase in tan f3. In the case where the evolution of mjj 
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i'm'jj^) is dominated by the term 3ffXt in their one-loop RGE [||,^ we see that 

- since increases with tan /? while ft is left essentially unaltered - the weak scale values 
of n^'jj^ and m'jj^ move closer to each other so that is reduced. As a result, for a larger 
tan P value the point will move out of the A-funnel region (modulo effects of the width of 
A) because niA becomes smaller than 2m^^. To return to the ^-funnel region, we must 
have a larger value of to also move rrij^^ to more negative values compensating for the 
reduction of niA with the increase in tan/?. For a fixed value of Aq, this means increasing 
mo, explaining why the Higgs funnel moves to the right as we increase tan/? in Fig. ^. A 
qualitatively similar behaviour can also be seen in Fig. but just in the last two frames 
since for fi > the Higgs funnel does not appear for the other choices of tan /?. 

A similar analysis can also help us to understand how the location of the Higgs funnel 
in the niQ — mi/2 plane depends on the choice of mt- For larger values of nit, and thus 
of ft, iTT-H^ evolves to more negative values while the evolution of fnjj^ remains essentially 
unaltered. Thus we move out of the ^4- funnel because now becomes larger than 2m^^, 
and to return to the A-funnel we must now reduce mo, so that the A- funnel (if it occurs) 
moves to smaller values of mo as mt is increased. Although we do not show figures here, we 
have verified that this is indeed the case for representative slices of the mo — mi/2 plane. 

Up to now in all our plots we have assumed j4o = 0. By changing the ^40 parameter, 
one is altering the intra- generation mixing between third generation sfermions, especially 
the top squarks. This mixing also reduces m^^. In Fig. we show the mo vs. mi/2 plane 
for tan/? = 10, Aq = —2 TeV and > 0. In this case, a forbidden region appears at low 
mo and mi/2) where the ti becomes the LSP. Along the edge of this region, a yellow/green 
band appears: the stop co-annihilation region, where the ti — Z\ mass gap is positive but 
quite small, so that Z\ can co-annihilate against t\ in the early universe, thus giving a relic 
density matching (1.1). We also see an /i-annihilation strip at low m;^/2 and mo ~ 2.75 TeV, 



where 2my ~ m/j and neutralino annihilations into SM fermions are resonantly enhanced. 



2.2 Models with scalar mass non- universality 

2.2.1 Generational non-universality: normal scalar mass hierarchy 

Motivation: The normal scalar mass hierarchy model (NMH) examines the effect of genera- 



tional non- universality in the SSB sfermion mass parameters 2C, 42]. While constraints 
from K — K mass difference restrict first and second generation scalar masses to be nearly 
universal, the constraints arising from B — B mixing are much less strict, and some non- 
universality of third generation matter scalars compared to first / second generation matter 
scalars can be allowed. In fact, it can be argued that the data actually favor such a 
case: we know that the measured value of BF{h — > S7) is in rather close accord with SM 
predictions, suggesting that the third generation sparticles that enter the 6 — > S7 loop 
diagrams are rather heavy- of order the TeV scale. Meanwhile, the 2 — 2>a discrepancy of 



^Here, Xt = ruQ^ + m?^ + mij^ + Al and Xb = rriQ^ + m?^ + m|-^ + Al. 

^This will be the case as long as squark masses and A-parameters are not simultaneously very small, 
since for the large values of tan/3 where the Higgs funnel occurs, the Yukawa couplings are typically larger 
than the electroweak gauge couplings. 
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mSDGRA : tanp-10, TeV. p>0, m^=171^ GeV 
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Figure 4: A plot of the toq vs. toi/2 plane in mSUGRA for Aq — —2 TeV, tan/3 — 10 with ^ > 
and mt — 171.4 GeV. The red-shaded regions are excluded by lack of correct EWSB or by presence 
of a charged LSP. Blue regions are excluded by direct SUSY searches at LEP2. Yellow and green 
shaded regions are WMAP-allowed, while white regions are excluded owing to ^^i^^ ^ 0.129. This 
plot includes a top-squark co-annihilation region adjacent to the excluded bulge at low mg and low 
mi/2 as well as an /i— annihilation strip at low mi/2 and mo ^ 2.75 TeV. Throughout the LEP2 
allowed region, rrih > 114 GeV. 

the measured {g — 2)^ against the SM prediction seems to favor rather light, sub- TeV scale 
smuon and muon sneutrino masses. A normal scalar mass hierarchy at the GUT scale with 
mo(l, 2) ^ mo(3) can reconcile the apparent tension between BF{b — > 57) and {g — 2)^ in 
SUSY models, and give a relic density in accord with (0). 

Parameter space: The parameter space of the NMH model is given by 

mo(l,2), mo, mi/2, ^0, tan/3, sign{fj,) (2.1) 

where mo(l,2) is the common GUT scale matter scalar mass parameter for first/second 
generation scalars at Mqut, while mo (3) = mH„ = it^h^ = "^o defines the remaining 
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scalar mass parameter, again at Q = Mqut- The Higgs scalar masses are taken here to be 
degenerate with mo (3), but could be independent as well. 

If we begin with a generic point in mSUGRA parameter space where fi^^/i^ » 0.129, 
and then dial mo(l,2) to successively lower values, the first/second generation slepton 
masses fall until they are low enough that bulk neutralino annihilation via light slep- 
tons and/or neutralino-slepton co-annihilation acts to reduce the relic density to WMAP- 
allowed levels. As a result of lowering mo(l,2), sleptons tend to be quite light. However, 
first/second generation squark masses are typically pulled up via RG running into the 
several hundred GeV to a TeV range. 

In Fig. ^, we show the ratio mo(l,2)/mo needed to reduce the relic density to the 
WMAP allowed value versus rriQ. We show results for three choices of mi/2 ■ 200, 300 
and 500 GeV. Solid curves are for tan /? = 10 while dashed curves are for tan f3 = 40. We 
fix j4o = and take fi > 0. At quite low mo, we are already in the stau co- annihilation 
region (bulk region for mi/2 = 200 GeV), so little or no reduction of mo (1,2) is needed. 
The curves terminate at the left because for still smaller values of mo, we hit the stau 
LSP region. As mo increases, a large reduction is needed to match the measured relic 
density, where ratios mo(l,2)/mo ~ 0.1 are common. At very large mo, no reduction in 
mo(l,2)/mo is again necessary as we enter the HB/FP region. 

2008/01/26 12.16 

NMH: fanp = 10& 40, Ag =0, \i >0, m, =171.4GeV 
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Figure 5: The value of mo(l, 2)/77io needed to bring various mSUGRA points into accord with 
the measured relic density versus mo, for Aq = 0, /i > 0, and tan (3—10 (solid lines) and 40 (dashed 
lines) in the NMH scenario. The three sets of curves correspond to toi/2 values of 200, 300 and 
500 GeV. 

Implications for collider searches: The very small first/second generation slepton 
masses in the NMH model imply that these will likely be directly accessible to LHC 
searches [^]. Even if this is not the case, branching ratios for chargino and neutralino 
decays to leptons via 2 or 3 body modes will be considerably enhanced, leading to SUSY 
cascade decay events at the LHC that are much richer in hard, isolated leptons than would 
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be expected in the mSUGRA model. Meanwhile, selection and smuon pair production - 
but not stau pair production - would likely be accessible to ILC searches. 

Implications for DM searches: While neutralino annihilation in the early universe is 
enhanced via light slepton exchange or slepton co-annihilation, squarks remain relatively 
heavy, and the neutralino is largely bino-like. Thus, both direct and indirect DM search 
predictions will be qualitatively similar to those generated in the mSUGRA model with 
n^oi'^) ~ "lo- For the case that we study in Table |l] below, we have explicitly checked that 
even the indirect detection signals at IceCube and Pamela detectors remain small despite 
the reduced masses of first/second generation sneutrinos and charged sleptons. 



2.2.2 Non-universal Higgs mass: one extra parameter case 

Motivation: In supersymmetric grand unified theories based upon the gauge group SO (10), 
the matter superfields of a single generation are contained in a 16-dimensional spinor rep- 
resentation of 50(10), ipiQ, which includes, in addition, a SM gauge singlet right-handed 
neutrino superfield. The Higgs superfields can be most simply accommodated in the funda- 
mental 10-dimensional representation <piQ. It is natural to expect that different multiplets 
would receive different soft masses at the GUT scale. Even if the soft masses for Higgs and 
matter scalars were common at some scale near Mp, RG running effects in the SO (10) 



theory would split the soft terms at Mqut (see Ref. |18] for explicit examples of soft term 
running in 50(10) SUSY GUTs). 

Parameter space: In the non-universal Higgs model with one additional parameter 
(NUHMl) the matter scalars receive a common squared mass parameter ttiq at Q = 
-^GUT) while both SU{2) Higgs scalar doublets Hu and acquire equal values for their 
SSB parameters, that are different from ttiq. Note that m|^^(= fnjj^) is just a parameter, 
not a physical mass squared, so its value can be either positive or negative (as can tuq [^] ). 
The parameter space is thus given by 

mo, 6^, mi/2, ^0, tan/3, sign{fi), (2.2) 

where = mo(l -|- S^) and = = sign{m^) ■ at the GUT scale. 

Given any parameter space point in the mSUGRA model with typically too high a 
relic density, one can always increase beyond its mSUGRA value of mo- A large value 
of > mo implies via the RGEs and EWSB minimization conditions a smaller weak 
scale value of and thus the possibility of mixed higgsino dark matter with a WMAP- 
allowed relic density (as in the mSUGRA HB/FP region), even though mo is not large. 
Alternatively, if is negative, the value of increases, but uia decreases: thus, by 
dialing m^ to a sufficiently negative value, we can get A-funnel annihilation, for any value 
of tan p. 

In Fig. ^ we show curves which illustrate the value of 6^ needed to move the mSUGRA 
relic density prediction into accord with (|l.lD , either by raising 6^ or equivalently, m^, as 
in frame a), or by lowering 6^ (and hence m,^), to negative values as in frame b). We show 
curves versus mo for mi/2 = 200, 300 and 500 GeV, and for Aq = 0, tan/3 = 10 and 40 
and ^ > 0. Curves terminate at both ends where we reach an end of the scanned mo space 
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or hit a forbidden region. At the low mo end, no diahng of is needed, since we are in the 
narrow stau co-annihilation region, or for mi/2 = 200 GeV, in the bulk region. When we 
move to larger mo, we eventually leave that region and large \5^\ values become necessary 
to lower 1^1 (or m^). As we continue to increase mo in the upper frame, smaller values of 
are necessary because of another assisting effect- the downward push of higgs mass-squared 
parameters from the top Yukawa coupling - is getting stronger. We smoothly reach the 
HB/FP region of mSUGRA where no dialing is required. The behaviour in the lower frame, 
where we adjust 6^ so as to hit the Higgs- funnel region is qualitatively different at larger 
values of mo: once we hit the HB/FP region, there is no need to have 6^ different from 
zero. The jump in the curves reflects the rapidity with which this region is reached. 
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Figure 6: The values of 5^ needed to bring various mSUGRA points into accord with the measured 
relic density versus mo for = and tan/3 = 10 (solid lines) and 40 (dashed lines) in the NUHMl 
model. The curves correspond to mi/2 values of 200, 300 and 500 GeV. 

Implications for collider searches: 

In the case with m^ > mo where we have MHDM (even though mo can be much 
lower than its typical HB/FP value in mSUGRA), the low value of implies that all the 
charginos and neutralinos will be quite light. Thus, they are more likely to be seen either via 
direct -ino pair production at the LHC, or to be produced at large rates in gluino and squark 
cascade decays. In general, for small \iJ,\, g and q cascade decay patterns become much more 
complex because many squark and gluino decay chains that are normally suppressed in the 
mSUGRA case become relevant. In addition, for small |//|, there is a smaller Z2 — Z\ mass 
gap, so Z2 2-body "spoiler decays" Z2 — > Z\Z and Z2 Zih are likely closed, and leptonic 
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3-body decays Z2 — > Z\li occur at observable levels. The dilepton mass edge from these 
three-body decays frequently serves as the starting point for sparticle mass reconstruction 
in SUSY cascade decays at the LHC |33]. Since is small, it is also possible that Z3 



also only decays via three-body channels. In this case, the dilepton mass distribution will 
contain three mass edges (though these may not all be observable), and its shape may 
provide further information about the nature of the neutralinos. 

In the case where < so that ~ 2m^^^ A as well as the other heavier Higgs 
bosons H and are much lighter than expected in mSUGRA, even at low-to-moderate 
tan /9 values. In this case - for a fixed value of tan (5 - direct detection of the heavier Higgs 
states ^4, H and at the LHC is more likely than in mSUGRA, and further, these states 
may also be produced in the gluino and squark decay chains via the decays of secondary 
chargino and neutralinos [45|. 



Implications for DM searches: 

For the case of ^ mg with small and MHDM, the enhanced higgsino component 
of the Zi leads to both enhanced direct and indirect DM detection rates compared to 
mSUGRA. This case has excellent detection prospects in the next generation of detectors 
such as XENON-100, LUX or mini-CLEAN. 

For the case where < with ~ 2m^^, the Zi remains nearly pure bino, so direct 
detection rates and indirect detection via neutralino annihilation to neutrinos in the solar 
core remain low, at values typical of mSUGRA models. However, indirect Zi detection 
via halo annihilations to gamma rays or anti-matter are all enhanced relative to mSUGRA 
(but not always to observable levels), since the halo neutralinos can still annihilate through 
the s-channel pseudoscalar resonance 



2.2.3 Non-universal Higgs mass: two extra parameters case 

Motivation: In SUSY GUT models based upon the gauge group SU (5), each of the MSSM 
Higgs superfields lives in different representations of the gauge group: G Ws, while 
G 7^5*. The two-extra-parameters NUHM model (NUHM2) [^] assumes independent 
Higgs field soft masses nn?^^ and f^'n^- The simplest assumption for the matter scalars is 
that they all acquire a common GUT scale mass mo, although they also would exist in 
separate 5* and 10 dimensional representations under SU{5). 

Parameter space: One form of parameter space for the NUHM2 model is 

mo, rriH^, mjj^, mi/2, Aq, tan/?, sign{fi) (NUHM2'). (2.3) 

However, the EWSB minimization conditions allow the new GUT scale parameters m'jj^ 
and m'jj^ to be traded for the weak scale parameters \fi\ and tua which are frequently easier 
to work with for phenomenological analyses: 

mo, ruA, mi/2, Aq, tan/3, sign{fi) (NUHM2). (2-4) 

Both forms of parameter space are allowed in Isajet spectra and event generation. In 
addition, Blazek et al. |47| adopt a Higgs SSB parameterization wherein the Higgs soft 
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masses are split evenly about a common mass which most generally is an independent 
parameter, miQ. Their parameterization is given by: 

rr^k,, = mlo (ItSh), (2.5) 

where 5h is dimensionless and can take either positive or negative values. If we choose 
mio = rriQ, we obtain a one-parameter extension that we refer to as the Higgs-splitting 
(HS) model. 

In Fig. 0, we illustrate the values of 6h needed to move the mSUGRA relic density 
prediction into accord with ( |1.1D by lowering 6h to negative values. This case gives rise to 
models with low /i and low uia in the HS model. We show curves of < needed to pull 



the mSUGRA relic density into accord with Eq. 1.1 versus mo for mx/2 = 200, 300 and 
500 GeV, and for tan/3 = 10 and 40, with = and fi > 0. Here the situation is similar 
to the NUHMl case in that less dialing is required from larger mo due to the increasing 
top Yukawa coupling effect. As in the previous figure, less dialing is needed for the larger 
value of tan /3. We also mention that if instead 5h is raised to large positive values (not 
shown in the figure), then instead one enters a WMAP-allowed region via ^l/P or ur/cr 
co-annihilation as discussed in Ref.[E2l. 



2007/10/29 14.05 




Figure 7: Values of Sh needed to bring various mSUGRA points into accord with the measured 
relic density versus toq for = and tan/3 — 10 (solid lines) and 40 (dashed lines) in the NUHM2 
HS scenario. The curves correspond to mi/2 values of 200, 300 and 500 GeV. 



Implications for collider searches: 

In the NUHM2 model, since /x and m^ are free parameters, one is free to choose both 
fj, and rriA to be small, so one can have MHDM and A-funnel annihilation contributions 
simultaneously. This type of model leads to the possibility of having light -inos and light 
A, H and at the same time. This would lead to a very complex and rich pattern of 
gluino and squark cascade decays at the LHC. 
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A new distinct possibility also arises in the NUHM2 model. In the MSSM scalar mass 
RGEs (see e.g. Ch. 9 of Ref. [^), the right-hand-side includes a term 

S = - + Tr [m^ - mi - Imfj + mjj + m|] 

which vanishes in the mSUGRA case, but is obviously non-zero if the Higgs soft masses 
are split. If the Higgs mass splitting is large, then the S'-term helps push the and 
soft terms to small values, leading to cases with left-slepton neutralino co- annihilation, or 
to cases with very light ur and cr squark masses, with squark-neutralino co-annihilation 
acting to reduce the relic density. The presence of very light left-sleptons or light ur, cr 
squarks, with other sfermions at the TeV scale, might be an indication of the HS scenario 
with large, positive 6h > 0. In this case, contrary to the prediction of many models, the 
lighter stau is dominantly fi- 

Implications for DM searches: 

For the NUHM2 model, all direct and indirect DM detection rates are strongly en- 
hanced if fj, is small, while only halo annihilation rates are enhanced if 2m^^ ~ mA- If 
instead parameters are in the region where agreement with the observed relic density occurs 
because the left-type sleptons are relatively light, we do not expect increases in either direct 
or indirect detection rates. Finally, in the light ilr, cr region, rates for direct detection 
of DM, as well as for its indirect detection via observation of high energy P^'s at IceCube 
will be enhanced since these depend mainly on neutralino-nucleon scattering via squark 
exchange, and ur squarks are light. Halo annihilation rates are also somewhat enhanced, 
since neutralinos can more easily annihilate into uu and cc pairs, giving rise to gamma ray 
and anti-matter signals [ p2[ . 



2.3 Models with non-universal gaugino masses 

Motivation: In mSUGRA, it is assumed that the gaugino mass parameters Mi, M2 and 
M3 unify to mi/2 at Q = Mqut- This holds true in supergravity models if the gauge 
kinetic function /ab ~ ^ABfihM), where A, B are gauge indices, and f{hM) is an arbitrary 
function of hidden sector fields /im; but common to all the gauge groups. More generally, 
the gauge kinetic function need only transform at the symmetric product of two adjoints. 
In this more general case, if the auxiliary field that breaks supersymmetry also breaks the 
grand unification gauge symmetry, GUT scale gaugino mass parameters need not unify |^8[ . 
Non-unified masses also occur in models of gaugino-mediated SUSY breaking and in 



various string-motivated models [50|. In models with mixed moduli-anomaly mediated 
SUSY breaking (MMAMSB) |Q, the gaugino masses are again split at Mqut, with the 
splitting proportional to the gauge group /3-functions. Motivated by these considerations, 
in the phenomenological models that we consider below, we will allow independent gaugino 
mass parameters at Q = Mqut- To isolate the effect, we will assume that just one of these 
mass parameters deviates from it unified value, and tune it to reproduce the measured 
value of the DM relic density leaving the other two at mi/2- 
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2.3.1 Mixed wino dark matter 

Parameter space: In mSUGRA, at the GUT scale one assumes Mi = M2 = ?tii/2 which 
leads to Mi ~ at the weak scale due to RG running. The Zi is usually a nearly pure 
bino state with too large a relic density. In anomaly-mediated SUSY breaking (AMSB) 
models |52|, the Zi is nearly pure wino-like, with too low a relic density. There exists, 



therefore, an intermediate situation with Mi M2 at the weak scale which gives the 
observed relic density, so that the Zi is mixed wino dark matter (MWDM) ||2^. Starting 
from the mSUGRA model boundary conditions at MguT; one can either increase Mi so 
that Ml > M2 = M3 = mi/2- 

mo, Ml, mi/2, Aq, tan/3, sign{fi) (the MWDMl case) (2.6) 

or lower M2 such that M2 < Mi = M3 = mi^2- 

mo, M2, mi/2, ^0, tan/?, sign{ij) (the MWDM2 case). (2.7) 

In the MWDMl case, since Mi, M2 and /i will all be much closer together, the Zi is 
actually a mixed bino-wino-higgsino state, while in MWDM2, since Mi ~ M2 <C /x, the Zi 
is more a pure bino-wino mixed state. 

In Fig. we show curves which illustrate the GUT scale ratio of ri = Mi/mi/2 needed 
to move the mSUGRA relic density prediction into accord with ( |1.1D . We show curves 
versus mo for mi/2 = 200, 300 and 500 GeV, and for tan/3 = 10 and 40. We take Aq = Q 
and /X > 0. At the lowest mo values ri = 1 since we are in the stau co-annihilation (or for 
mi/2 = 200 GeV in the bulk) region. For mo values beyond this, we have to dial Mi so as to 
obtain a mixed bino-wino-higgsino Zi to be in concordance with the observed relic density 
measurement. For very large values of mo approaching the HB/FP region, the required 
value of ri once again begins to reduce as long as m^^ > Mw so that ZiZi — > W~^W~ is 
accessible in the early universe . The up-turn in the mi/2 = 200 GeV curves at large values 
of mo occurs when this reaction becomes kinematically suppressed: in this case, a larger 
value of ri once again allows this reaction without which the relic density tends to be too 
large. However, for mi/2 = 200 GeV and tan/3 = 10, and mo > 1.45 TeV, we are deep 
enough into the HB/FP region so that even for mSUGRA the annihilation cross section 



via Z* exchange is large enough to get agreement with (1J-) even with m^^ < Mi^/.^ In 
the case of the curves for mi/2 = 300 and 500 GeV, m^^ always remains above M^y so that 
the HB/FP region of mSUGRA is smoothly reached. The dashed curves exhibit analogous 
behaviour. 

Implications for collider searches: 

In MWDM models with Mi ~ M2 such that the relic density constraint ( |1.1[ ) is fulfilled, 
the mass gap m^^ — m^^ generally tends to be ~ 20 — 30 GeV for MWDM2 models (larger 
for MWDMl, where the neutralino also has a higgsino component), somewhat smaller than 
the mass gap of ~ 50 GeV found in models with MHDM (such as the HB/FP region of 
mSUGRA). This means that at collider experiments, the same- flavor/opposite-sign isolated 



*This portion of the curve is excluded by the LEP constraint on mfr; 



Wi ■ 
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Figure 8: Values of ri — Mi/mi/2 needed to bring various mSUGRA points into accord with the 
measured relic density versus mo for = and tan/3 = 10 (solid lines) and 40 (dashed lines) in 
the MWDMl scenario. The curves correspond to mi/2 values of 200, 300 and 500 GeV. For the 200 
(300) GeV solid curves, the regions with mp 1-4 (1.94) TeV are excluded because the chargino 
is too light. 

dilepton mass spectrum should have a single visible edge around 20-30 GeV arising from 
Z2 — > ZiU decays. Moreover, the shape of the dilepton spectrum should correspond to one 
where the neutralino eigenvalues have the same sign. This should be distinct from MHDM 
models which tend to have the higher mass edge and possibly a shape corresponding to 
opposite signs of the neutralino mass eigenvalues |^3[. Also since MHDM models have a 
small /i parameter, edges due to decay may also be visible. The small Z2 — Z\ mass 
gap suppresses 3-body decays of Z2 more than 2-body decays: then the branching fraction 
for the radiative decay Z2 Zij, which is normally very suppressed in mSUGRA, can 
reach the 10% level in the MWDM models Finally, when Mi is raised, it feeds 

into raising masses of especially right-type sleptons relative to their mSUGRA predictions. 
Likewise, when M2 is reduced, left- squark and slepton masses get reduced relative to 
mSUGRA. 

Implications for DM searches: 

In the MWDMl scenario, the Zi becomes a mixed bino-wino-higgsino state, and as 
a result all direct and indirect DM detection rates are boosted relative to mSUGRA - 
sometimes by an order of magnitude or more due to the enhanced higgsino component. 
In the MWDM2 scenario, where the Zi develops a smaller higgsino component, direct 
detection and indirect detection rates are only slightly enhanced, while indirect DM 
detection rates from halo annihilations can again be boosted by an order of magnitude or 
more, but not necessarily to observable levels. 

2.3.2 Bino-wino co-annihilation (BWCA) 

Parameter space: In the BWCA scenario where Mi ~ —M2 at the weak scale, there 
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is very little mixing between the bino and neutral wino states even when these are very 
close in mass. If |Mi| is just slightly smaller than IM2I, and |//| is relatively large, the 
lightest neutralino remains bino-like, but because m^^ ~ , bino- wino co-annihilation 

processes in the early universe reduce the relic density to the observed level. In the BWCA 
case, the parameter space is the same as in the MWDM case, except that Mi and M2 are 
now opposite in sign. As with the MWDM case, one can either raise \Mi\ by about a factor 
2, or lower IM2I by about the same factor to attain —Mi ~ M2 at the weak scale. 

In Fig. ^ we show curves which illustrate the GUT scale ratio r2 = M2/mi/2 needed 



to move the mSUGRA relic density prediction into accord with (1.1). We show results 
versus mo, for mi/2 = 200, 300 and 500 GeV, and for tan/? = 10 and 40. We take ^0 = 
and fi > (upper frame) and ^ < 0, the sign favored by the E821 (g — 2)^ experiment 
(lower frame), at least for lower ranges of mo and mi/2- That the ratio r2 remains close 
to —1/2 in the upper frame, for all mo values in between the stau-coannihilation and the 
HB/FP regions, reflects the fact that the evolution of gaugino masses does not depend on 
sfermion masses at 1-loop. The solid lines for tan/3 = 10 in the lower frame show very 
similar behaviour, except that /U is now negative. In contrast, there is a large flat region at 
low values of mo in the tan /3 = 40 case shown by the dashed lines. We have traced this to 
the fact that the Higgs funnel region has already opened up even for tan/? = 40, and that 
the funnel region is contiguous to the stau co-annihilation (and for mx/2 = 200 GeV, also 
the bulk) region. 

Implications for collider searches: In the BWCA scenario, the Z2 — Z\ mass gap 
becomes very small: of order 15-30 GeV typically. Thus, as in the MWDM case, one 
expects a m{fJ) mass edge to be visible owing to Z2 ZiU decay, since all Z2 two-body 
spoiler decay modes are kinematically closed. The mass gap should be much smaller than 
that typically expected from mSUGRA models. The m{i^i~) distribution for opposite- 
sign/ same flavor dileptons should contain a single mass-edge, "one hump bump", unlike 
the case of MHDM, which favors a "two-hump-bump" since Z3 would be light as well.^ 

The small mass gap (required for effective co-annihilation) strongly favors two-body 
decays over three-body decays and results in large branching fraction for the radiative 
decay Z2 — > ^17: reaching over 30% in the scans presented in Ref. [^]. If Z2 were at 
rest, the 7 from the radiative decay would be mono-energetic but rather soft. However, for 
fast moving Z2 secondaries from gluino and squark cascade decays, the photon energy gets 
boosted, and in the BWCA case there should be an observable signal also in the multi-jet 
plus isolated photon plus E!^^^^ channel at the LHC: see Fig. 16 of Ref. [p^]. 

At e^e~ colliders, the relative minus sign between Mi and M2 leads to enhanced 
production of .^^1^2 pairs compared to mSUGRA predictions, and also the predictions in 
the MWDM case. This provides a way of distinguishing between the BWCA and MWDM 
frameworks which otherwise have a very similar mass spectrum. Moreover, operating just 
above .^^1^2 threshold, one might make detailed studies of Z2 decay branching fractions, 
including the radiative mode Z2 — > Zij, which will result in events with almost mono- 
energetic single photons recoiling against "nothing". 



^We thank Dr. Theodore Geisel for coining related expressions. 



- 21 - 



2007/10/25 10.15 



BWCA2: fanp = 10& 40, Ag =0, n >0, =171.4GeV 




200 



300 



500] 



0.2 0.4 0.6 0.8 



1.2 1.4 1.6 1.8 2 
m„(GeV) 



2008/01/21 18.25 

BWCA2: fanp = 10& 40, Ag =0, ii <0, m, =171.4GeV 



g~ -0.5 

§ -0.6 
II 

^ -0.7 
-0.8 
-0.9 
-1 




500 



0.2 0.4 0.6 0.8 



1.2 1.4 1.6 1.8 2 
mg(GeV) 



Figure 9: Values of r2 — M2/mi/2 needed to bring various mSUGRA points into accord with 
the measured relic density versus mo for = and tan/3 =10 (solid lines) and 40 (dashed lines) 
in the BWCA scenario. The upper frame is for /i > while in the lower frame we take /i < 0. The 
curves correspond to mi/2 values of 200, 300 and 500 GeV. 



Implications for DM searches: In the BWCA case, since the Zi remains purely bino- 
like, rates for indirect DM detection remain low, similar to results from the corresponding 
mSUGRA case in many instances. Rates for direct DM detection can be far below the 
sensitivity of any proposed detector if the sign of Mifj, is negative, i.e. for fi > in the 
BWCA2 case, and for ^ < in the BWCAl case, because of cancellations in neutralino 



couplings that enter the direct DM detection rate calculations |24]. 



2.3.3 Low IM3I dark matter: compressed SUSY 

Parameter space: The low IM3I dark matter (LM3DM) scenario arises by starting with 
mSUGRA parameter space, but (for mo < 1 — 2 TeV) lowering the GUT scale value of 
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jMal relative to Mi = M2 = mi/2 |55, 25|. Lowering IM3I results in smaller gluino and, via 
RGE effects, also squark masses. These effects feed into the MSSM RGEs and affect the 
running of fn^jj^ - effectively diminishing the downward push from the top quark Yukawa 
coupling - resulting in lower \^\ values, and hence MHDM. The MHDM case can be easily 
compatible with the observed relic density constraint since there is enhanced neutralino 
annihilation to WW, ZZ and Zh states in the early universe. Thus, the parameter space 
is given by 

mo, 'mi/2, M3, Aq, tan/3, sign(^) (the LM3DM case). (2-8) 

Here, M3 can be of either sign. Although the first and second generation sfermion masses 
are essentially unaffected by the sign flip, m^^, and through /i, also chargino and neutralino 
masses, do show clear dependence on the relative sign between M3 and Mi^2- 

In Fig. H^, we show the GUT scale ratio of = M3 / ^2 needed to move the mSUGRA 



relic density prediction into accord with (LI) versus mg. As always, we show results for 



mi/2 = 200, 300 and 500 GeV, and for tan/3 = 10 and 40. We take Aq = and fi > 0. 
For the solid curves, we see that once we are away from the stau co-annihilation (and 
for the mi^2 = 200 GeV case, also bulk) region, we need to reduce |M3(GUT)| to obtain 
the correct relic density. Since in this scenario we are lowering //, the degree of dialing is 
generally smaller for larger mo due to increasing top Yukawa coupling effects, just as in the 
NUHMl model with positive 5(f,. The situation is more complicated for the dashed curves 
where tan/^ = 40. For M3 < 0, the Higgs funnel already starts to appear for this relatively 
low value of tan/3. For the red and blue dashed curves corresponding to mi/2 = 200 and 
300 GeV, respectively, this funnel region is contiguous with the bulk/stau co-annihilation 
region, so that remains at -1 for uiq < 600 GeV. For yet larger values of mo, |r3| needs 
to be dialed down though, because of the proximity of the Higgs funnel, by not quite as 
much as for the corresponding tan/3 = 10 case, until the HB/FP region is reached. For 
the Tni/2 = 500 GeV curve, the Higgs-funnel region occurs for 450 GeV < mo < 800 GeV, 
and is well separated from the stau co-annihilation region. Thus some tuning of (but 
again, not as much as in the tan /3 = 10 case) is needed for mo values away from the very 
narrow stau co-annihilation region, and again for large mo values outside the Higgs-funnel 
region, until the HB/FP region is reached at mo ~ 1.9 TeV. For M3 > 0, although the 
Higgs-funnel does not occur for tan/3 = 40, s-channel A/H exchange does significantly 
enhance neutralino annihilation amplitudes for mo ^ 2mi/2- ^ result, the value of 
needed varies more slowly at the low mo end for the dashed curves than for the solid curves. 



We should mention that a related scenario, dubbed "compressed SUSY", has been 
suggested by Martin ||5^. In compressed SUSY, M3 is lowered, but also by choosing large 
negative values of the Aq parameter, a rather light ti state can be generated. Then if 
m^^ > nit, neutralino annihilation via ZiZi — > tt, which does not suffer the usual P- 
wave suppression because of the large top quark mass, can dominate in the early universe, 
resulting in the observed relic density in a different way. The phenomenology of compressed 
SUSY models has been examined in Ref. [^7| . 

Implications for collider searches: Since nig and rriq are lowered relative to m^ values. 
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Figure 10: Values of the GUT scale ratio ra = 17111/2 needed to bring various mSUGRA 
points into accord with the measured relic density versus toq for Aq = and tan/3 = 10 (solid lines) 
and 40 (dashed lines) in the LM3DM scenario. The curves correspond to mi/2 values of 200, 300 
and 500 GeV. 



the Tevatron and LHC reach is enhanced compared to corresponding mSUGRA predictions. 
In mSUGRA, the chargino mass limit rn^^ > 103.5 GeV from LEP2 usually pre-empts 
the limit from direct gluino searches at the Tevatron because in models with gaugino mass 
unification, rrig ~ 3.5m^^, and thus nig > 350 GeV, which is not far from the reach of a 
2 TeV pp collider. However, in LM3DM with non-unified gaugino masses, the gluino can 
have a mass as low as ~ 200 GeV, while the chargino remains in the LEP2 allowed region, 
m^^^ > 103.5 GeV . Thus, a significant chunk of LM3DM parameter space is open to gluino 
pair and gluino-squark searches at the Tevatron |^8| . Since the gluino mass is lowered with 
respect to m^^, the radiative gluino decay g — > Zig may be the dominant decay mode of 
the g, with a branching fraction as high as 85% in regions of parameter space where 
squarks are very heavy. In this case, the decay of a gluino leads to a single high pT jet, 
and gluino pair production will look more like squark production. 

Since the Zi is MHDM in the LM3DM scenario, the Z2 — Z\ mass gap is again lowered 
compared to mSUGRA predictions, and the Z2 ZiU decay should be visible at LHC 
searches. In the LM3DM case, the mass gap is typically of order 30-80 GeV, and there 
should be contributions to the m{ii) distribution as well, thus allowing the MHDM 
scenario to be distinguished from BWCA or MWDM. The large higgsino content of the 
neutralino Zi will also tend to lead to a large branching ratio for gluino decays to third 
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generation quarks so that the reach of the LHC will be enhanced by tagging 6-jets 

For ILC searches, it is more likely that squark pair production (especially titi produc- 
tion) will be accessible. In addition, since the Z3, Z4 and W2 states are lighter than the 
corresponding mSUGRA parameter space points, it is more likely that various -ino pair 
production reactions would be accessible to ILC searches, allowing the reconstruction of 
chargino and neutralino mass matrices. 

Implications for DM searches: In the LM3DM scenario, since we expect both light 
squarks and low n (leading to MHDM), direct DM search rates are enhanced relative to 
mSUGRA by up to two orders of magnitudel Rates for events at IceCube are enhanced 
by up to three orders of magnitude! Similar enhancements are seen in gamma ray and 
anti-matter search predictions arising from neutralino annihilation in the galactic halo. 
Thus, the LM3DM scenario seems a boon for direct and indirect DM searches. 

2.3.4 High IM2I dark matter: left-right split SUSY 

Parameter space: In the high IM2I dark matter scenario (IIM2DM), the parameter set is 
the same as in the MWDM2 scenario, except that now the SU (2) gaugino mass parameter 
M2 is dialed up in magnitude ||2^. This increased value of IM2I feeds into the MSSM RGEs 
by first pushing rrij^^ to higher values (than in the universal gaugino mass case) during its 
evolution from Q = Mgut- Then as RG running continues, the top Yukawa coupling terms 
take over, and m'\j^{Q) begins to be reduced. Since its positive peak- value was higher than 
in the canonical case with universal gaugino masses, m?^^ attains a relatively small negative 
value when the evolution is stopped at the weak scale. Finally, since //^ ~ — m|^^(weak), 
we end up with a small weak scale parameter, and MHDM. 

In the HM2DM scenario, the large value of |M2(Mgut)|) via RG evolution, lifts the 
SSB masses of SU{2) doublet matter scalars to large values, so that left-sleptons, and to 
a smaller extent also left-squarks, are much heavier than right- ones. Thus, the IIM2DM 
model can be regarded as left-right split SUSY. In this model, all light third generation 
matter sfermions {irif^) then tend to be predominantly right- states, whereas in most 
models, the 61 tends to be mainly 6^, and ti tends to be a mixed left-right squark state. 
Further, since |/i| is small, the HM2DM model leads to a spectrum with light Zi, Z2, 
and Wi states, which tend to be higgsino-like or mixed bino-higgsino. The W2 and are 
nearly pure winos, and also very heavy. 

In Fig. 11, we show values of the GUT scale ratio of r2 = M2/mii2 needed to bring the 



neutralino relic density prediction into accord with (LI). We show curves versus mo for 
mi/2 = 200, 300 and 500 GeV, for tan/3 = 10 and 40, with ^0 = 0. Since Aa^^^^ oc ^iM2 
we take sign{fj,) = sign{M2) so their product is positive in accordance with the measured 
value of {g — 2)^ pO| . As in the LM3DM model, flipping the sign of M2 (without changing 
its magnitude) causes the A-funnel to open up for tan f3 = 40. But now the funnel region 
remains merged with the stau-coannihilation region even for mx/2 = 500 GeV, so that for 
the lower range of mo no additional dialing is necessary and r2 remains at —1. 

Implications for collider searches: In the IIM2DM scenario, since we again have 
MHDM, the Z2 — Z\ mass gap tends to be in the range 30-80 GeV, so that the mass 
edge from Z2 Z\U decays, and possibly also from Zj, decays, should be seen in gluino 
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Figure 11: Values of the GUT scale ratio r2 = M2/TO1/2 needed to bring various mSUGRA 
points into accord with the measured relic density versus toq for Aq = and tan/3 = 10 (solid lines) 
and 40 (dashed lines) in the HM2DM scenario. The curves correspond to toi/2 values of 200, 300 
and 500 GeV. 



and squark cascade decay events at the LHC. The shape of the distributions may also make 
it possible to glean information about the relative signs of the neutralino mass eigenvalues. 
Also, the increased higgsino content of the lighter states should again lead to increased 
6-jet multiplicity in SUSY events at the LHC. 

At the ILC, the production of chargino and neutralino pairs would vary in a contrasting 
way compared to mSUGRA because the low lying gaugino states would be essentially devoid 
of wino components. If third generation squarks and sleptons are accessible to ILC searches, 
then a variation in beam polarization would reveal all these states to be predominantly 
right-type states, and their pair production cross sections would decrease with increasingly 
left-polarized beams. 

Implications for DM searches: Since in the HM2DM scenario, the Z\ is a mixed bino- 
higgsino state, its signal for spin- independent direct detection should be observable at the 
next generation of detectors (super-CDMS or 100 kg noble liquid detectors) over much 
of the parameter space. The signal at IceCube or Antares would also be boosted 
by up to two orders of magnitude compared to mSUGRA and may be observable over a 
substantial portion of parameter space. Rates for detection of gamma rays and anti-matter 
from neutralino halo annihilation are also boosted relative to mSUGRA by 1-2 orders of 
magnitude. 
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3. Illustrative benchmark cases 



In this section, we list some benchmark cases of models with universality and non-universality. 
We start with the mSUGRA model, and adopt a point in parameter space with 

mo, mi/2, Ao, tan/?, sign^fi) = 300 GeV, 300 GeV, 0, 10, +1 

with rrit = 171.4 GeV. This point is listed in column 2 of Table |^. We see that we get 
il^^/i^ = 1.1 for this point which is conclusively in conflict with (p..l|), and so excluded 
assuming standard Big Bang cosmology and thermal relic neutralinos. For every other 
model in Table ||, we relax the universality assumption and allow one additional parameter 
that we tune to bring the model into the DM-allowed range with ^^_^h'^ ~ 0.1. The 
point here is to be able to compare and contrast the spectra along with other features 
of each of these DM-allowed models with the corresponding spectrum of the mSUGRA 
model, and with one another. We also list at the bottom the Isatools output of h? , 
BF{b 57), the SUSY contribution to the muon anomalous magnetic moment Aaj^^'^^, 
the spin- independent cross section for the elastic neutralino-proton scattering asi{Zip), 

and the Higgsino content of the neutralino Rfj = \J + ^2"*^'*^ in the notation of Ref . Q . 
We use IS A JET 7.76 to generate this Table. 

In the first non- universal case, the NMH model, we see that the first/second generation 
SSB parameter mo(l,2) = 54 GeV in order to obtain the observed relic density. Then, 
the and states (and also the JIr and JIl states) have much reduced masses than the 
corresponding mSUGRA case, while nif^ ^ > ^"^^ ^.Iso gluino, chargino and neutralino masses 
are essentially unchanged from their mSUGRA values. The low value of mg^ = 128.9 GeV 
and 10.5 GeV mass gap between cr/ hr and Zi ensure a high rate for neutralino annihilation 
and co-annihilation in the early universe. In addition, while BF{b — >■ sj) remains near the 
measured and SM value because third generation squarks and charged Higgs bosons are 
heavy, the value of Aa^^^^ is enhanced, thus reconciling these two possibly disparate 
measurements. 

The second non-universal case, labelled NUHMl^, comes from the NUHMl model 
where is dialed up to 549 GeV so that /i becomes small and we have mixed higgsino 
DM, even though mo is far smaller than in the mSUGRA HB/FP region. In this case, 
has been reduced so much that the point is actually LEP2 excluded. Note that Rjj 
has risen to 0.84, signaling a higgsino-like Zi. The branching ratio for the decay b ^ 
is slightly reduced. The salient feature is the direct detection cross section, which is now 



in the range of well-tempered neutralino models [pSj, 31 1 with MHDM, is 37 times higher 
than the corresponding mSUGRA value, and very close to the current 90% CL limit from 
Xenon- 10 search |61]. 



The third non- universal case, labelled NUHMl^, dials to —728 GeV, which raises fx 
to large values but lowers niA to be just above 2m^^ so that yl-funnel annihilation reduces 
the relic abundance, even though tan /3 is not large. In this case, sfermions and lighter -inos 
have essentially the same masses as in mSUGRA, but now we have relatively light Higgs 
boson states A, H and accessible to LHC searches. Also, BF{b — > 57) is somewhat 
enhanced due to the light entering the tH^ loop contribution to this decay. 
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The fourth non-universal case is the HS model and has the two Higgs soft masses split 
about the common value mg, and has both MHDM and some ^-funnel annihilation, with a 
light spectrum of Higgs bosons, as well as charginos and neutralinos. All the Higgs bosons 
and all the charginos and neutralinos should be accessible at an electron-positron collider 
operating at a center of mass energy just above 500 GeV. The higgsino component of the 
Zi, which is larger than in mSUGRA, leads to considerable enhancement of the direct 
detection cross section. 



Model 


mSUGRA 


NMH 


NUHMl,, 


NUHMIa 


HS 


parameter 




mnfl, 2) 




' ' "(p 


5h 

^ ±1 


special value 




54 


549 


-728 


-1.36 




385.1 


386.5 


105.8 


748.5 


269.3 




729.7 


722.1 


731.4 


733.4 


728.9 




720.8 


658.4 


724.3 


720.5 


720.1 




523.4 


526.5 


484.1 


624.5 


505.8 




656.8 


659.8 


642.2 


689.5 


645.4 




364.5 


216.2 


364.8 


365.8 


373.4 


rrieji 


322.3 


128.9 


322.5 


321.9 


301.8 


mfi 


317.1 


317.6 


317.8 


316.4 


299.3 


vv 2, 


411.7 


412.7 


264.7 


754.8 


321.1 




220.7 


219.5 


91.1 


234.9 


196.6 




412.5 


413.5 


268.1 


754.6 


322.9 




391.3 


392.7 


137.3 


747.1 


277.1 




220.6 


219.4 


117.4 


234.5 


198.1 


^Z, 


119.2 


118.4 


69.0 


121.5 


115.4 


ruA 


520.3 


521.9 


584.5 


268.5 


279.0 


mH+ 


529.8 


531.4 


593.8 


281.6 


292.0 


rrih 


110.1 


110.1 


109.8 


110.5 


109.8 


BF{b 57) X 10"^ 


1.1 


0.10 


0.11 


0.11 


0.10 


3.0 


3.1 


2.5 


4.3 


3.4 


Aa^ X lO^o 


12.1 


27.2 


17.9 


9.3 


13.7 


asi{Zip) X 109 (pb) 


2.1 


2.1 


78 


1.2 


27 




0.15 


0.14 


0.84 


0.06 


0.26 



Table 1: A comparison of the characteristics of mSUGRA with corresponding characteristics in 
models with scalar mass non-universality that lead to the observed relic abundance of DM. Input 
parameters and resultant sparticle masses in GeV units, together with the predicted neutralino 
relic density, BF(b —^ sj) and Aa^, the SUSY contribution to the anomalous magnetic moment 
of the muon, the direct detection cross section for the Zi, and finally, the higgsino content of the 
Zi. In each case, we fix mg = 'mi/2 — 300 GeV, Ag = 0, tan/3 = 10 and rrit — 171.4 GeV, and for 
each non-universal model tunc the parameter in the first row to its special value shown in row 2 to 
reproduce the observed relic abundance. 

In Table |, we continue this comparison with the same mSUGRA point in Table ||, but 
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this time for models with non-universal gaugino mass parameters. We show the results for 
this mSUGRA point once again for the convenience of the reader. In column 3, we consider 
the MWDMl model where we raise the GUT scale value of Mi to 490 GeV, resulting in 
a Z\ state that is a mixed bino-wino-higgsino state. The heightened wino and higgsino 
components of Z\ allow for enhanced Z\Z\ — > W^W~ in the early universe, thus putting 
the model into accord with the measured DM abundance. The model has direct detection 
rates enhanced by a factor ~ 7 over mSUGRA, to the 10~^ pb range. Except for Z\ which 
now has a mass of 195 GeV compared to 119 GeV in mSUGRA, the sparticle spectra are 
almost the same in the two cases. The heavier Z\ state implies that the Zi — Z\ mass gap 
is about 29 GeV compared to ^ 100 GeV in mSUGRA, and that its higgsino-content is 
somewhat larger than in mSUGRA. 

In the next column we consider a BWCA scenario where we dial Mi(Mgut) to 
—480 GeV to obtain the observed relic density. Again, we see that except for m^^ which 
moves to 202 GeV, sparticle and Higgs boson masses as well as the higgsino content of 
the Z\ are essentially the same as for the corresponding mSUGRA model. The W\ — Z\ 
mass gap is now just 18 GeV, so bino-wino co-annihilation acts to reduce the relic density, 
even though the Z\ remains in a nearly pure bino-like state. The small Z^ — Z\ mass gap 
compared to the value of m^, which might be deduced at LHC, would signal a model with 
non- universal gaugino masses. Notice that for reasons detailed in Ref. |31], the Z\ direct 
detection cross section is far lower than any proposed detector can probe. 

In column 5, we attain the observed relic density by dialing M3(Mgut) from 300 GeV 
to 160 GeV to obtain a viable LM3DM model. We see that we now have a much lighter 
spectrum of sparticles than in mSUGRA, not only squarks and gluinos, but also charginos 
and neutralinos (but not sleptons). In this case, we would expect huge sparticle production 
cross sections at LHC and complicated cascade decay chains. Since the reduction of M3 
hardly affects slepton masses, sleptons are not much lighter than squarks even though the 
sfermion mass scale is just 300-400 GeV. The Z\ has a significant higgsino content leading 
to MHDM with ?i, Z^ — Z\ mass gap of 58 GeV, while the — Z\ mass gap is only slightly 
below M^, so 3-body Z3 decays will also occur, though the dilepton mass will be peaked 
close to Mz- The large higgsino content also results in a correspondingly enhanced Z\ 
direct detection cross section - 7 x 10~^ pb - close to the 90% CL limit from Xenon-10 
DM searches. We also see that the relatively light top squarks and charginos significantly 
reduce BF{h — > S7). 

Finally, in the last column in Table we show results for the HM2DM model where we 
get agreement with (LI) by raising M2(Mgut) to 900 GeV. We see that this gives a very 
low value of \jl so that the Z\ becomes MHDM with Rfi = 0.67. We see that the left-type 
squarks and especially the left-type sleptons are now much heavier than in mSUGRA, while 
qPi and Ir are hardly affected, leading to the left-right split spectrum referred to earlier. 
The lighter chargino, and the three lightest neutralinos are all lighter than in mSUGRA, 
while is very heavy. The Z2 — Z\ mass gap is ~ 47 GeV, while the Z3 — Z\ mass 
gap is just 64 GeV. The branching ratio for 6 ^ 57 is reduced from its mSUGRA value 
because t\ and the higgsino-like chargino are relatively light, while the DM direct detection 
cross-section is large, as is typical of models with MHDM. 
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15 
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72 


34 




0.15 


0.25 


0.16 


0.50 
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Table 2: A comparison of the characteristics of mSUGRA with corresponding characteristics in 
models with gaugino mass non- universality that lead to the observed relic abundance of DM. Input 
parameters and resultant sparticle masses in GcV imits, together with the predicted neutralino 
relic density, BF{b sj) and Aa^, the SUSY contribution to the anomalous magnetic moment 
of the muon, the direct detection cross section for the Zi , and finally, the higgsino content of the 
Zi. In each case, we fix mo = mi/2 = 300 GeV, Aq = 0, tan/J = 10 and mt = 171.4 GeV, and for 
each non-universal model tunc the parameter in the first row to its special value shown in row 2 to 
reproduce the observed relic abundance. 

4. General characteristics of relic-density-consistent models 

In this section, we abstract general features of the various models that we have introduced 
earlier by performing scans over model parameters, where we keep only parameter points 
which lead to a relic density ^^i^"^ ~ 0.11. Wc also reject models with < 103.5 GeV 
from LEP2 searches. In the case of the mSUGRA model, we scan parameters over the 
range: mo : - 5 TeV, mi/2 : - 2 TeV, Aq = 0, tan/9 = 10,30,45,50,52,55 and both 
signs of 11. Thus, our scans will include the stau co-annihilation region, the HB/FP region 
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at large rriQ and the ^-funnel at large tan/3. For models with non-universality, in order 
to have manageable parameter space scans, we restrict ourselves to a scan over mSUGRA 
parameters m-o : — 2 TeV, mi/2 : — 1.5 TeV, = 0) tan/3 = 10 and fj, > (except 
for the HM2DM and BWCA2 models with negative M2 where we take /U < since this is 
somewhat favored by the measurement of {g — 2)^). For models with non- universality, for 
every set of mSUGRA parameters in our scan, we adjust the special additional parameter 
listed in the first rows of Tables |^ and |2| to bring the neutralino relic density h'^ into 
accord with (|1.1| ). Our upper limit on is chosen somewhat arbitrarily to avoid too 

much fine-tuning. We allow much larger values of mo just in the mSUGRA region, since 
in this case it has been argued that fine-tuning in the HB/FP region is not very large |ll|| . 
We limit ourselves to lower values of ttiq < 2 TeV in the models with non-universality since 
MHDM characteristic of the low values in the HB/FP region can be attained for all 
values of uiq. 

4.1 Implications for collider searches 

In Fig. |l^ we show the value of rrig vs. rriu^ (as a representative value of the approximately 
degenerate squark mass) for both signs of mSUGRA, as well as for eight other models 
with non-universal SSB terms indicated in the legend on the figure. Each dot shows the 
gluino and up-squark mass for a model, with parameters chosen so that the neutralino 
relic density saturates (|l.l| ). Notice that in this figure, as in several subsequent ones, not 
all the colors are visible since some model points are overwritten by other model points. 
The diagonal dashed line for m^^ = nig shows that when we require that the observed 
relic density be obtained, with the exception of the two branches in the HB/FP region of 
mSUGRA where niQ (and hence the squark mass) is very large, all models yield mg ~ rriq. 
The two dotted lines denote the approximate reach of the CERN LHC with 100 fb^^ of 



integrated luminosity if rrig ~ rriq, as adapted from Ref. |62]. Most of the models scanned 
lie within reach of the CERN LHC, which is partly an artifact from the upper limit we 
take on mi/2 in our parameter space scans. For the two branches in the HB/FP region of 
mSUGRA shown by the red and blue dots at large tuq in the figure, experiments at the 
LHC will be sensitive to models where rrig ^1-8 TeV. 



In Fig. [Ij, we show these DM-allowed models in the m^^ vs. uih plane. Here, we note 
a clear trend in all models: heavier ti squarks are correlated with larger values of rrih, 
largely because top- Yukawa radiative corrections to nih increase with the stop mass. For 
many models with ^ Mz, then h ~ i^sM so that the LEP2 lower bound of 114.1 GeV 
would be applicable. We have not required this bound in our analysis for reasons discussed 
earlier, and also to be able to show the trend of m/j with other observables. The lightest 
value of rrij^ occurs in the LM3DM model, which allows m^^ as low as ~ 200 GeV, although 
even lighter values of m^^ (readily accessible even at the Tevatron |6^) would be allowed 
if we also admitted variation of ^o- 

In Fig. we show the lighter chargino mass m^^^^ vs. — m^^. We denote by 
the dashed line the region where — m^^ < Mz. Below this line, the spoiler decay 
modes Z2 — > ZiZ or Zih are kinematically closed, so that Z2 must decay via 3-body 
modes like Z2 — > Ziii. The mass edge in the invariant dilepton mass distribution from this 
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Figure 12: Predictions for rrig vs. rriu^ from various models witli = 0, mt = 171.4 GeV 
and ^ > (except in the cases with the blue dots for the mSUGRA model, light blue dots for the 
BWCA2 model and magenta dots for the HM2DM model with M2 < for which we have < 0), 
but where the special parameter in the various non-universal mass models has been dialed to yield 
ri^^/i^ ~ 0.11. We fix tan/3 — 10 except for the mSUGRA model where we allow tan/3 = 10, 30, 
45, 50, 52 and 55. The approximate 100 fb~^ reach of CERN LHC is denoted by the dotted lines, 
while a dashed line denotes where m^^ = m,g. Here, and in subsequent figures, dots for some of 
the models are covered up by other dots, and so are not visible. 



decay p3[ |, which can serve as the starting point for sparticle mass reconstruction at the 
LHC as discussed earlier, will be visible as long as its branching fraction is not strongly 
suppressed |3^. From the figure, we see that most of the models that give rise to the 
correct relic density also predict that the spoiler modes are closed so that the m{U) mass 
edge will hkely be visible! Exceptions arise from the stau-co-annihilation and yl-funnel 
regions of mSUGRA, or at vestiges of other models which already had the correct relic 
density (because the mSUGRA parameters were in the stau co-annihilation region) so that 
no special non-universality parameters needed to be adjusted to get the correct relic density 
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Figure 13: Predictions for mh vs. ■mi_^ from various models with Aq ~ Q, mt — 171.4 GeV and 
the sign of as in Fig. ^ but where the special parameter of non-universal mass models has been 
dialed to yield Slg ~ 0.11. We fix tan/3 = 10 except for the mSUGRA model where we allow 
tan/3 = 10, 30, 45,' 50, 52 and 55. 



{e.g. the gray points from the MWDMl model). Notice also that there are models where 
the mass gap is very small. Even for these models it is likely that the Z2 will be sufficiently 
boosted in its production via decays of much heavier squarks/gluinos so that the daughter 
leptons have large enough transverse momenta so as to be detectable. 



In Fig. 15, we show predictions for rriyj^ and m^j for WMAP-allowed models. The 
approximate reach of the ILC500 {^/s = 500 GeV) and ILCIOOO (with ^ = 1000 GeV) 
are shown by the dashed and dotted lines , respectively, which delineate the kinematic 
limit for W^W^ or f^f^ pair production. Here, we see that it is quite easy to evade the 
ILC reach and still have a neutralino relic density consistent with ( |l.lD . This is in contrast 
to prejudices from studies in the mid-1990s which favored the bulk annihilation region of 
mSUGRA, which then implied sparticle mass ought to be quite light, and likely accessible 
to LEP2 and ILC500 searches B. The upper bands of mSUGRA model parameter points 
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Figure 14: Predictions for vs. m^^ ^''^Zi ^oi^^ various models with Aq — Q,mt — 171.4 GeV 
and the sign of /i as in Fig. |l2|, but where the special parameter of non-universal mass models has 
been dialed to yield — 0.11. We fix tan/3 = 10 except for the mSUGRA model where we 

allow tan/3 = 10, 30, 45, 50, 52 and 55. The dashed line denotes the point where m^^ ~ 
where the two-body spoiler decay mode Z2 ^ Z\Z turns on. 

correspond to the HB/FP region, while the bands of points at low but high fn^^ 
correspond to the stau co-annihilation region in mSUGRA or in MWDMl models. 

4.2 Implications for (5 — 2)^ and BF(h — > 57) 

The rare decay b ^ s'j has always been interesting for SUSY (as well as other new physics) 
studies, because the SM and the new physics contributions both occur at the one-loop 
order, and so are likely to be comparable if the particles in the new-physics loop have 
masses of about the weak scale. This is indeed the case for weak scale SUSY. The branching 
fraction BF[b — > 57) has been measured by the CLEO, Belle and BABAR collaborations; a 
combined analysis finds the branching fraction to be BF{h — > 57) = (3.55ib0.26) x 10~^, 
while a recent SM prediction [^] finds BF{h — > 57) = (3. 15 it 0.23) x 10~^. The theoretical 
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Figure 15: Predictions for vs. m^^ from various models with Aq = 0, rrit = 171.4 GeV and 
the sign of as in Fig. ^ but where the special parameter of non-universal mass models has been 
dialed to yield fl^ — 0.11. We fix tan/3 = 10 except for the mSUGRA model where we allow 
tan P = 10, 30, 45, 50, 52 and 55. The dashed lines denote the approximate reach of ILC500, while 
dotted lines mark the approximate reach of ILCIOOO. 



error in the SUSY case may be somewhat larger. In Figure 16, we show predictions for 
BF{b — > 57) in SUSY models where ~ 0.11, against the value of nig. We see 

that for models with low nig, large deviations from the SM prediction are likely, although 
cases in agreement can be readily found. As nig increases, the SUSY loop contributions to 
BF{b S7) are suppressed. In the absence of an underlying theory of flavor, we should 
be careful in drawing strong inferences from this figure since even a small amount of flavor- 
violation in the textures of SSB parameters could significantly alter these predictions, 
with little impact on implications for direct searches at the LHC. That BF[b — > 57) has 
potentially larger SUSY contributions for mSUGRA models than in non-universal mass 
models is, of course, an artifact of our scans: for the mSUGRA model we scan large values 
of tan f3 while we fix tan /3 = 10 for models with non-universal masses. 
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Figure 16: Predictions for BF{h sj) vs. rrig from various models with Ao = 0, nit = 171.4 GeV 
and the sign of /i as in Fig. |l2|, but where the special parameter of non-universal mass models has 
been dialed to yield ^^_^h^ — 0.11. We fix tan/3 = 10 except for the mSUGRA model where we 
allow tan/3 = 10, 30, 45, 50, 52 and 55. The dotted line denotes the central value of the combined 
experimental measurements, while the dashed line denotes the corresponding SM prediction. The 
SUSY contribution to the branching fraction is sensitive to tan (3 which is varied for the scans of the 
mSUGRA model, but fixed at tan/3 = 10 for the scans in the case of non-universal mass models. 



Recent measurements of the muon anomalous magnetic moment show an apparent 
deviation from SM predictions. Combining QED, electroweak, hadronic (using e^e~ — > 
hadrons to evaluate hadronic loop contributions) and light-by-light contributions, and com- 
paring against measurements from E821 at BNL, a positive deviation in = ^'^^ of 

Aa^ = al'^P - ^ = 22(10) x 10-^° (4.1) 

is reported in the Particle Data Book [40|, i.e. a 2.2a effect.^ 

®More recent analyses|^^ report a larger discrepancy if only electron-positron data are used for the 
evaluation of the hadronic vacuum polarization contribution; the significance of the discrepancy is, however, 
reduced if tau decay data are used for this purpose. 
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One-loop diagrams with Wi — and Zi — /ii^2 in the loop would give super symmetric 
contributions to a^, perhaps accounting for the (rather weak, yet persistent) discrepancy 
with the SM. In Fig. |l^, we show Aa^^^"*^ versus m^^ . The dashed line indicates the central 
value of the experiment/theory discrepancy as presented by the Particle Data Group. We 
see that a variety of models are able to account for the discrepancy as long as rrij^^ < 1 — 1.5 
TeV for mSUGRA, but only about 500 GeV in the case of models with non-universal mass 
parameters. This is partly a consequence of the fact that in mSUGRA our scans include 
large values of tan (3 in mSUGRA but are limited to tan /3 = 10 for non-universal models, 
bmce 

Aa^uSY ^ tan/?, had we allowed larger values of tan (3 in our scans of models with 
non-universality, then the Aa^ projections would increase beyond those plotted here, and 
consistency with the present central value would be possible for values of second generation 
slepton masses beyond the reach of a 1 TeV ILC. 

4.3 Implications for direct detection of dark matter 

In Fig. [l^, we show the spin-independent neutralino-proton scattering cross section, cal- 
culated with IsaReS program from the IsaTools package, versus fn^^. A significant 



uncertainty in the cross section comes from the value of the pion-nucleon S-term |68D . In 
this plot we assumed S = 45 MeV, but larger values can increase our predictions by factor 
of about three. This plot is an update of similar results presented in Ref. ^l| in that 
it includes additional models. We also show the current limit established by the Xenon- 
10 collaboration ||6l| (solid line), along with the projected reaches for the SuperCDMS 



(25 kg) [||] (dashed line), LUX 300 kg [^] (dot-dashed line) and Xenon-1 ton |71[ (dot- 
ted line) experiments. The reach contours have been generated assuming a standard local 
density and velocity profile. 

We see two distinct classes of models. In the first class, the neutralino-nucleon cross sec- 
tion falls off with , while in the second class - models with a well-tempered neutralino 
with significant higgsino component - this cross section asymptotes to about 10^^ pb, 
within the reach of the next generation of detectors such as LUX-300 kg, Xenon-100 or 
super-CDMS. It is important to realize that this second class includes several of the specific 



models that we have considered where agreement with (1J-) is obtained via a significant 
higgsino component in the Zi so that we have either mixed higgsino DM or mixed wino- 
bino-higgsino DM. This higgsino component then leads to a large cross section for Zip 
scattering via diagrams involving h and/or H exchanges, where the Higgs bosons couple 
to the proton via both its quark and its gluon content. The neutralino annihilation rate in 
the early universe generally falls off with increasing m^^, so that for heavier neutralinos, a 
larger higgsino content is necessary to maintain the relic density at its observed value: it is 
precisely this increased higgsino-content that maintains the direct-detection cross section 
around 10^^ pb even for large values of in the upper branch of the figure. There are, 
however, many models where accord with the observed CDM relic density is obtained by 
adjusting the masses to get either stau co-annihilation or bino-wino co-annihilation or Higgs 
funnel annihilation. In these cases of the bino-like LSP, the direct detection cross section 
falls with m^^ to below the sensitivity of even It noble element detectors for neutralino 
masses below about 400 GeV. There even are cases with rrii? < 100 GeV where - due to 
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Figure 17: Predictions for Aaf^^^^ vs. m^^ from various models witli Aq = 0, nit ~ 171.4 GeV 
and the sign of /i as in Fig. |l^, but where the special parameter of non-universal mass models 
has been dialed to yield H.^ ~ 0.11. We fix tan/3 = 10 except for the mSUGRA model where 
we allow tan/3 = 10, 30, 45, 50, 52 and 55. The dashed line denotes the central value of the 
measured deviation from SM expectations as reported by the Particle Data Group, though some 
recent analyses would infer an even larger deviation as discussed in the text. Note that Aa^ is 
sensitive to tan/3 which is varied for mSUGRA, but for scans of the non-universal mass models, is 
fixed to be 10. 



interference between various contributing amplitudes |72, ^] {e.g. mSUGRA with fi < 0) 
- the neutralino-nucleon scattering cross section drops to well below 10^^*^ pb, which is 
below the projected sensitivity of all proposed detectors to date. These same interference 
effects frequently lead to a reduced cross section in the BWCA2 case where we also take 
/X < 0, the sign favored by the value of Aa^. 

The neutralino may also scatter inside a detector via its spin-dependent coupling to 
the nucleon due to its couplings to the Z or to squarks. In Fig. 19 we show how this 
spin-dependent cross section is expected to scale with the corresponding spin-independent 
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Figure 18: Predictions for asi{Zip) vs. m^^, generally regarded as the figure of merit for direct 
detection experiments, in various models with Aq — 0, mt — 171.4 GeV and the sign of fj, as in 
Fig. but where the special parameter of non-universal mass models has been dialed to yield 
^^^h? ~ 0.11. We fix tan/3 = 10 except for the mSUGRA model where we allow tan/3 = 10, 30, 
45, 50, 52 and 55. We also show the projected reach of selected direct detection experiments. 



cross section in the various models that we have considered. It is striking to see that 
while the spin-dependent cross section in relic-density-consistent models may be as low 
as 10~^ pb, in well-tempered neutralino models where agreement with the relic density is 
obtained by adjusting the higgsino content of Zi this cross section is always larger than 
10~^ pb well above the projected sensitivity, (Tsd(^ip) ^ 4 x 10~^ pb, of the proposed 
COUPP experiment with a target mass of It |73|. This is of course, because higgsinos have 
a large coupling to the Z-boson. We note that there may be an observable signal via spin- 
dependent couplings even for cases where the prospects for direct detection via the spin- 
independent neutralino interaction appear to be hopeless. We also remark that the 50 kg 
prototype of the COUPP detector is projected to have a sensitivity asoiZip) ~ 4x 10~^ pb. 
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Spin -de pendent vs. Spin -independent: Cross Section 




Figure 19: Predictions for <jsd{Zip) vs. asi{Zip) in various models with Aq — Q,mt — 171.4 GeV 
and the sign of /i as in Fig. |l2|, but where the special parameter of non-universal mass models has 
been dialed to yield — 0.11. We fix tan/3 — 10 except for the mSUGRA model where we 

allow tan/3 = 10, 30, 45,' 50, 52 and 55. 



4.4 Implications for indirect detection of dark matter: neutrino telescopes 

In Fig. we show the flux of muons with > 50 GeV which is expected from neutralino 
capture by the sun, with subsequent neutralino annihilation in the solar core to states ( 
for some recent work, see Ref. |75]). To calculate these and subsequent indirect dark matter 
detection rates, we use the DarkSUSY |76] - Isajet interface. The flux of high energy 
from neutralino annihilation in the solar core depends on both the neutralino capture cross 
section as well as on the neutralino annihilation cross section. The capture rate mainly 
depends on the spin- dependent neutralino-nucleon cross section, the main contribution to 
which comes from Z-exchange processes that are enhanced if the neutralino has a significant 
higgsino content. Note that it is easily possible to have a large spin-independent neutralino- 
nucleon scattering cross section, and yet a small spin-dependent cross section, leading to 
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undetectable rates in the IceCube experiment This is exemplified in Fig. 2C, we again 
find that the well-tempered neutralino models aggregate around an asymptotic regime of 
~ 10 — 100 events/km^/yr. The approximate reach of IceCube for > 50 GeV is around 
the 40 events/km^ /yr [^^, so again many of the models with mixed gaugino/higgsino dark 
matter stand a good chance of indirect detection via neutrino telescopes, whereas models 
where (|1.1[) is satisfied in other ways fall below the detectable level. We remark though 



that if ~ 100 GeV, even the signal from BWCA2 may be in the detectable range. We 
mention that projections for IceCube, unlike those for indirect detection from neutralino 
annihilation to anti-matter or gamma rays (discussed next), are only slightly sensitive to 
the DM halo profile. 




Figure 20: Predictions for muon flux with > 50 GeV $^ vs. m^^ from various models with 
Aq — 0, nit — 171.4 GeV and the sign of fi as in Fig. 12, but where the special parameter of 
non-universal mass models has been dialed to yield ^zj^^ — O-H- We fix tan/? = 10 except for the 
mSUGRA model where we allow tan/3 = 10, 30, 45, 50, 52 and 55. The region above the dashed 
line denotes the approximate reach of the IceCube neutrino telescope. 
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4.5 Implications for indirect detection of dark matter from halo annihilations 



An alternative method for indirect detection of dark matter is to search for the debris 
resulting from dark matter annihilation in the galactic halo. One promising method is 
to search for GeV-scale 7-rays, which could come directly from Z\Z\ 77 via a box 
and triangle diagrams, or via Z\Z\ qq, where hadronization and decay lead to gamma 
rays via the g — 7r° ^ 7 chain. The direct process occurs at low rates, but would have 
a characteristic signal at the source with ~ ''^Zi^ whereas gamma rays coming from 
quark hadronization should be more abundant, but will yield a continuum distribution in 
with a cut-off at ni^^ . 

Since 7-rays from neutralino annihilation should propagate undeflected through the 
galaxy, a good place to look is the galactic center, where the DM density is expected to 
be high. In Fig. we show the flux of 7-rays coming from the direction of the galactic 
center with Ej > 1 GeV, in units of events/cm^/s. The result is very sensitive to the 
choice of the galactic dark matter density profile, as well as to (unknown) details of how 
clumpy the halo distribution is. We show results for the Adiabatically Contracted N03 halo 
profile |7^, where the deepening of gravitational potential wells caused by baryon in-fall 
leads to a higher DM concentration in the center of the Milky Way, and a concomitantly 
larger gamma ray flux. Other halo distributions, such as the Burkert profile (7^], where 
the central DM halo cusp is smoothed out by significant re-heating, predict gamma ray 
fiuxes that may be four orders of magnitude smaller] The reach of the GLAST satellite 
experiment is indicated by the dashed line The important point is that models with 
large s-wave neutralino annihilation cross sections cluster around an asymptotic level of 
~ 10~'^/cm^/s, while models which rely on co-annihilation such as BWCA or MWDM2 
predict much lower gamma ray fluxes. A notable difference between signals from halo 
annihilation versus signals from direct and neutrino detection is that the halo annihilation 
signals can be enhanced by moving 2m^^ onto the ^-resonance if neutralinos have 
enhanced annihilation through the A- funnel in the early universe, then they can also readily 
annihilate through the ^-funnel in the galactic halo (this does not hold true for the h 
and H resonances, which are dominantly p-wave, or the Z pole, which is not resonance 
enhanced [^). In the case of 7-ray signals, we see the orange dots from NUHMl/^ model 
now populate higher rate levels than the BWCA and MWDM2 cases, whereas for direct 
and signals in Figs. and 2C, the NUHMl^ signal was comparable to or even lower 
than the BWCA and MWDM2 models. 

Another characteristic signature of DM halo annihilations is the detection of large 
fiuxes of anti-particles such as ps, e^s or anti-deuterons D. For positrons and antiprotons, 
we evaluate the averaged differential antiparticle fiux in a projected energy bin centered at 
a kinetic energy of 20 GeV, where we expect an optimal statistics and signal-to-background 
ratio at space-borne antiparticle detectors |81, 82]. We take the experimental sensitivity 
to be that of the Pamela experiment after three years of data-taking as our benchmark. 
For Ds, we evaluate the average differential anti-deuteron fiux in the 0.1 < Tg < 0.25 GeV 
range, where Tg stands for the anti-deuteron kinetic energy per nucleon, and compare it to 
the estimated GAPS sensitivity for an ultra- long duration balloon-borne experiment |83| 
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Figure 21: Predictions for gamma ray flux with > 1 GeV from various models witli Aq — 0, 
mt = 171.4 GeV and the sign of /i as in Fig. |l^, but where the special parameter of non-universal 
mass models has been dialed to yield fi^ ~ 0.11. We fix tan/3 = 10 except for the mSUGRA 
model where we allow tan/? = 10, 30, 45, 50, 52 and 55. We adopt the Adiabatically Contracted 
N03 DM halo profile. The region above the dashed line denotes the approximate reach of the 
GLAST experiment. The flux predicted from less cusped halo proflles may be down by as much as 
four orders of magnitude. 



(see Ref. |8J] for an updated discussion of the role of antideuteron searches in DM indirect 
detection) . 

In Fig. ^ we show the flux of ps assuming the Adiabatically Contracted N03 halo 
profile; results from using the Burkert profile yield results typically a factor of 10-20 below 
these. The models with mixed higgsino dark matter cluster at high levels of around 
10~^ events/GeV/cm^/s/sr while the A-funnel annihilation case of NUHMl^ populates the 
10~^ — 10~^ events/GeV/cm^/s/sr range. The co-annihilation cases BWCA and MWDM2 
and stau-co-annihilation in mSUGRA lie at much lower levels. 



In Fig. 23 we show the flux of e"''s using the Adiabatically Contracted N03 halo profile; 
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Figure 22: Predictions for anti-proton flux from various models with Aq = 0, mt — 171.4 GeV 
and the sign of /i as in Fig. |l2|, but where the special parameter of non-universal mass models has 
been dialed to yield fi^ ~ 0.11. We fix tan /? = 10 except for the mSUGRA model where we allow 
tan/3 = 10, 30, 45, 50, 52 and 55. We adopt the Adiabatically Contracted N03 DM halo profile. 
The region above the dashed line denotes the approximate reach of the PAMELA experiment. 

results from using the Burkert profile yield results about a factor of 3-5 lower^. Our projec- 
tions are not optimistic. The models with mixed higgsino dark matter cluster at the ~ 10~^ 
events/GeV/cm^/s/sr level, which may be just below the Pamela reach |^5[. The A-funnel 
annihilation case of NUHMIa populates the 10"^*^ — 10~^ events/GeV/cm^/s/sr range. 
The co-annihilation cases BWCA and MWDM2 and stau-coannihilation in mSUGRA are 
again at much lower levels. 



In Fig. 24, we show the predicted flux of anti-deuterons expected in a kinetic energy 



^To reach earth before losing too much energy and annihilating, the positrons must originate from 
annihilation much closer to earth than for ps or 7s; thus, predictions for their flux are less sensitive than 
those for anti-protons and gamma rays to the choice of halo profile. Different halo distributions mainly 
differ on the DM density near the galactic center, but agree on the local DM density. 
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Positron Detection : Ad. Contr. N03 HM 
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Figure 23: Predictions for positron flux from various models with Aq — 0, mt — 171.4 GeV and 
the sign of as in Fig. ^ but where the special parameter of non-universal mass models has been 
dialed to yield Slg — 0.11. We fix tan/3 = 10 except for the mSUGRA model where we allow 
tan/3 = 10, 30, 45, 50, 52 and 55. We adopt the Adiabatically Contracted N03 DM halo profile. 
The region above the dashed line denotes the approximate reach of the PAMELA experiment. 

range Tg = 0.1 — 0.25 GeV using the Adiabatically Contracted N03 halo profile, suitable for 
detection by the proposed GAPS experiment. Results using the Burkert profile tend to be a 
factor of 10-20 lower, about the same as for anti-protons. Models with mixed higgsino dark 
matter populate the 10~^^ events/GeV/cm^/s/sr level, and should be accessible to GAPS 
via the long duration balloon flight. The NUHMl^i model populates points just below to 
just above the GAPS sensitivity level, while the co-annihilation models give results which 
are generally beyond reach of any foreseeable experiment. 

In drawing inferences for prospects for indirect detection from the scans of non- 
universal mass models discussed in this section, we should keep in mind that we have 
shown results for just the Adiabatically Contracted N03 halo profile (for detailed compar- 
ison of halo profiles, see e.g. Ref. [^] and 1^^) and fixed tan/3 = 10. Typically, we have 
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Anti-deuteron Detection : Ad. Contr. N03 HM 




Figure 24: Predictions for anti-deuteron flux from various models with Aq = 0, mf = 171.4 GeV 
and the sign of ji as in Fig. |l^, but where the special parameter of non-universal models has been 
dialed to yield ^^i^'^ — O-H- We fix tan/3 — 10 except for the mSUGRA model where we allow 
tan/3 = 10, 30, 45, 50, 52 and 55. We adopt the Adiabatically Contracted N03 DM halo profile. 
The region above the dashed line denotes the approximate reach of the GAPS experiment. 

found that the indirect searches are most sensitive to the higgsino component in the Zi. 
It is important to note that direct detection and indirect detection via halo annihilation 
both grow as tan/3 is increased. 

5. Summary and concluding remarks 

If the observed cold dark matter Q is interpreted as thermal relic neutralinos of i?-parity 
conserving supersymmetric models, then the determination of the CDM relic density ( |l . l| ) 
provides a very strong constraint, effectively reducing the dimension of model parameter 
space by one unit. It is then reasonable to ask how this relic density measurement constrains 
on what other experiments searching for SUSY might or might not observe. 
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Indeed many groups have analyzed the imphcations of the measured value of the CDM 
relic density for SUSY signals at the LHC within the mSUGRA framework. Toward the 
end of Sec. ||, we enumerated several broad conclusions that were drawn from these studies. 
In order to test the robustness of these conclusions, it is necessary to examine how these 
are affected if we relax the untested universality assumptions that are the hallmark of the 
mSUGRA framework. Motivated by this, as well as by the fact that most of the relic- 
density-allowed range of parameters lies on the periphery of mSUGRA parameter space, 
we have examined a variety of models where universality of scalar or gaugino SSB mass 
parameters is relaxed via the introduction of just one additional parameter that is then 
adjusted so that the thermal relic density of neutralinos matches (|L1|), by adjusting either 
the neutralino composition or its mass. In Sec. ^, we show explicit examples of these various 
models that lead to broadly similar sparticle mass spectra, and compare and contrast the 
features of the different models with the paradigm mSUGRA framework, and with one 
another. 

Prior to the analyses of non- universal models, there were several prejudices inferred 
from studies based on mSUGRA, and frequently held to be true, including: 

1. The relic-density-consistent bulk region implies a variety of light sparticles accessible 
at the LHC, and possibly the ILC; 

2. The Higgs-funnel region only occurs at large tan /3, where down-type Yukawa cou- 
plings are necessarily large, so that sparticle decay cascades are modified, with con- 
comitant effects on collider signatures; 

3. The higgsino-content of the neutralino LSP can only be large enough to get agreement 
with only if scalars are essentially decoupled at the LHC; 

4. The lighter 6-squark is dominantly bL while the lighter stau is dominantly tr. 

We have seen that even in relatively innocuous one-parameter extensions of mSUGRA each 
of these conclusions is false. For instance, the HS model allows rapid neutralino annihilation 
via light ur/cr with other sparticles much heavier, or via fi which is dominantly tl, the 
Higgs-funnel occurs for any value of tan (3 in the NUHMl model, and we can have MHDM 
for rather small scalar masses, also in the NUHMl model. While it is definitely worthwhile 
to correlate the implications of one observation with what might and might not be seen 
in other experiments, our analysis highlights the fact that such inferences are frequently 
dependent on underlying assumptions. In particular, we caution against drawing broad 
conclusions about what is or is not likely at the LHC based upon studies of just the 
mSUGRA model. 

In Sec. we have performed scans over the parameter space of the mSUGRA as well 
as over eight of its one-parameter extensions to abstract features common to relic-density- 
consistent models. We end by summarizing our broad conclusions based on these scans. 

• In mSUGRA, a well-tempered neutralino LSP can only be obtained in the HB/FP 
region, where squark and slepton masses are far heavier than the lightest charginos. 
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neutralinos and gluino. In non-universal models, we can easily have a well-tempered 
neutralino with niq ~ rrig. Indeed except for the HB/FP region of mSUGRA, squark 
and gluino masses are typically comparable in relic-density-consistent models. In a 
similar vein, we also note that while Higgs-funnel enhancement is possible only for 
very large values of tan /3 in the mSUGRA framework, if we allow for non-universality 
of Higgs SSB parameters, we can have the Higgs funnel for any value of tan /3. 

In many relic-density-consistent models, the Z2 — Z1 mass gap is usually less than Mz, 
so that two-body spoiler decays modes of Z2 are kinematically closed. This means 
that at least one dilepton mass edge (and perhaps more) is likely to be detectable 
at LHC. The location of the dilepton mass edge(s) is a rather clean signature of 
supersymmetric models, and often serves as the starting point for sparticle mass 
reconstruction. 

Most relic-density-consistent models should lead to observable signals at the LHC. 
In contrast, while models where accord with the observed relic density is obtained by 
tempering the higgsino-content of the neutralino will likely be accessible at a 1 TeV 
electron-positron collider, in other scenarios sparticles may simply be too heavy to 
be accessible. 

In well-tempered neutralino models, the mechanism that enhances annihilation in 
the early universe also tends to enhance the direct DM detection rate. In particular, 
models tempered via the higgsino content of the LSP typically have a si (Zip) ~ 10~^ 
pb, which ought to be accessible to the next set of direct detection experiments, in- 
cluding LUX, Xenon-100, WARP, mini-CLEAN and SuperCDMS: see Fig. These 
experiments may also provide a measure of the mass of the halo DM particle(s), as- 



suming that it is not very heavy compared to the target nucleus [86].^ If a signal is 
found in these direct detection experiments, it can be directly compared to expec- 
tations based on SUSY model parameters extracted in experiments at the LHC and 
especially the ILC to test whether thermally produced neutralinos indeed saturate 
the measured value of the cold DM density [^, or whether DM, like visible matter, 
turns out to have more than one component. 

• Likewise, these models have elevated rates for indirect DM detection via neutrino 
telescopes. In this case, the flux of muon neutrinos tends to be above ~ 10 
events/km^ /year for E^, > 50 GeV. In many such models, the signal should be acces- 
sible at the IceCube detector, as can be seen from Fig. 

• Finally, well-tempered neutralino models also have elevated rates for indirect DM 
searches via neutralino annihilation in the galactic halo into gamma rays and anti- 



matter, especially if the higgsino component is enhanced: see Fig. 21-24. These rates 
have large uncertainties associated with the presently unknown galactic dark matter 
density profile. But if a signal is found, it can be compared to expectations using 



^Direct detection experiments with different target nuclei ranging over a wide range of masses may thus 
provide clear evidence for multiple WIMP components in the galactic halo. 
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model parameters extracted from LHC and ILC measurements, and the measured 
halo annihilation rate can be used to determine the DM halo profile |87|. 
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